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INTRODUCTION 


There has been a recognized need for realistic steady and unsteady transonic 
airload prediction methods suitable for flutter, gust load, and static aero-' 
elastic analysis of flight structures (particularly wings and control surfaces) 
since the first advent of high speed, compressibility-induced control problems 
of World War II aircraft (References 1-3). Although the level flight speeds 
at that time rarely exceeded 400 mph, the contemporary thick, unswept wings 
occasionally encountered serious supercritical flow problems in high speed 
dives (Reference 1, Chapter 9). 

The advent of the turboprop and turbojet engine led to an entirely new 
generation of fighters, bombers, ASW aircraft, and both military and civilian 
transport aircraft operating at high subsonic Mach numbers. Transport air- 1 
craft, in particular, are characterized by thick, moderately swept wings of 
large aspect ratio. These characteristics are required for efficient, long 
range cruise. Also, adequate high-lift performance is needed to provide 
efficient airport performance within noise limits and with minimum wing 
structural weight to obtain a maximum ratio of useful load to gross weight. 

The first generation of jet transport aircraft, in fact, have been operating- 
with supercritical wing flow at the high speed cruise point. Typically, one 
finds local steady flow Mach numbers, referred to the swept chord, ranging 
from 1.2 to as high as 1.6 at high angles of attack. This is the typical 
situation for a 35° swept wing, for example, at a flight Mach number of . 85 , 
where the freestream Mach number component normal to the quarter chord line 
is only .7. Thus, one often sees a two-fold increase of local Mach number 
in the "supercritical", (locally supersonic) flow region. 

The current generation of turbofan-powered wide-body jets operate in essen- 
tially the same speed and supercritical flow regime. Primarily, the new 
turbofan powered aircraft provide improved fuel economy and important 
reductions in airport jet noise, rather than increased speed. It is 



- 1 predicted that future technology advances in steady flow aerodynamics may| 

i be used to allow: (l) increases of wing thickness, (2) reductions of sweep, 

I ■ 1 

(3) increases of cruise speeds, or (4l) STOL (powered high lift) capability. 

I 

--! i 

In view of the current stress on economy and noise reduction, the first two 

— i 

possibilities may gain greater emphasis than speed increases for CTOL, and 

especially for STOL aircraft; however 1 , such trade-offs are a proper function 

| ( 

i of mission analysis. Suffice it to say that the aircraft evolved from any 

; of the. above described advances in technology, including STOL aircraft, will 

I yield designs whose wings will operate well into the .supercritical flow ! 

1 i 

regime. Furthermore, military aircraft must be flutter certified to I 

I V = 1.15 V- (V is design or structural limit speed), and commercial aircraft 
D D I 

, must certify to 1.2 V^, which further extends the transonic domain of | 
interest to the aeroelastician. In military aircraft the current emphasis 

! i 

is upon transonic maneuvering at high C T , to allow higher load factors under 

I ■* J ] 1 

; buffet-free conditions. j 1 

1 _ I __ 

For the purposes of the present study', the extensive body of experimental! 

I 1 

data and recent calculation procedures for steady flow can be regarded as ( 

available for defining a non-uniform steady flow field "environment" (with 

\ | 1 

variable thermodynamic state variables) into which small amplitude, unsteady! 

, I | 

aerodynamic flow perturbations must propagate. The aeroelastically induced 

perturbations are small compared to the steady flow variables; therefore, J 

it is' appropriate to linearize the equations of motion with respect to the 

: i 

unsteady perturbations. Since the steady flow field varies, this approach . 

I i 

is often referred to as "local linearization" in the literature (References ; 
4-7)- While this is an appropriate! description of the mathematical 1 

process, the term "local linearization" is sometimes identified with methods 

: I 

employing additional simplifying assumptions with respect to the boundary 1 J 

conditions, especially neglect of sonic line reflections in the super- ^ j 

critical region. The removal of this' deficiency, which is prevalent in all 1 

I . 

of the current literature on unsteady transonic flow, will be a primary , 
feature -of -the -present analysis.— ______ — 1 — | 


l 
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Also, incorporated in the present method is a series of suggestions for ! 
possible simplifying key assumptions | for idealizing the steady flow, based 
on extensive steady flow research within the aerospace industry with 
regard to transonic airfoils and wings, including , the results of airfoill 
research, and experience derived from aerodynamic wing design and develop- 
ment of transport aircraft. • 

I 

The above mentioned key assumptions will facilitate the application of 1 

I 1 

other previously well developed pieces of technology in: (l) unsteady aero- 
dynamics for purely subsonic and purely supersonic flow, and (2) relevant 
technology in acoustic propagation through non-uniform (layered) moving I 
media. 


In summary, the analytical approach outlined below and in various Appendices 
is believed to be one which is physically well grounded, and which will 
appeal to practical flutter analysts |because it will: 

o' Autpm¥t’iT:aiLly _ TFduCe~to'" s'iandard results “under subcriffca'l 


flow or purely supersonic 


flow conditions. 


o Should cause only a nominal increase in computation time when j 

efficiently implemented on the digital computer. I 

I 1 

o Will require only a few readily available input parameters to l 

characterize the essential features of the supercritical steady 

flow environment. 

i 

o Is capable of continuous refinement, as improved- steady flow I 
field data description becomes more routinely available. 
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SYMBOLS 


speed of sound 

transfer matrix, Eq (E-4);; also coefficients in Eqs (E~9^>E-95) 
wing reference length (semichord), Eq (D-3) 1 

= \JyP~-~L) Fig Al, Eqs (A-35j C-3) ; also, transfer matrix ,Eq (E-4)' 


coefficient in Eqs (E-94, 


E-95) 


chord (Fig 3)j section coefficient (lift, drag or pitching moment 
(Fig 3) | > 

coefficients in velocity potential solution, Eqs (A-15, E-2) 1 

differential operator (d/dz) 1 

I 1 

base of natural logarithm I 

I I 


— _ transfer- f unct -i-on-or-c oe ffie-i-e nt-,— w-i t h- di-ff erent-meani-ngsy- 


-!- “4“ 


depending on special subscripts or supercripts. See Eqs (A-20 

to A-32 , e-24, e-26, E- 30I31) I 


frequency 


Hz 


particular integral forcing functions, ,Eq (E-5) 


Green's function transfer 


matrix for shock induced downwash 


) 


-* -f 


A 


at wing surface, Eqs (B-20, B-21); also Green's function transfer 

| I 

matrix for shock impingement disturbances due to W , Eqs B-27 |to 

A j 

B-29; also other meanings] depending on subscripts and superscripts, 
Eqs (B-38, B-42, B-46) I 

complex number operator | 

1 

unit vectors ; j 

integral, Eq (C-19) j 

inverse of matrix A, Eq (E-25) j 


* > 



Dl 3 


unit vector (see above); also reduced frequency, Eq D-3 

'| > 

aerodynamic downwash induction matrix^ Eqs (B-3, B-4-9); also^ 
wave parameter, Eqs (A-10, A-ll, 0 ,- 2 .) 

Laplace transform operator Eq (E-79) 

Mach number 
unit normal vector 

lifting pressure; absolute static pressure (Fig 2) 


interference downwash induction matrices ,Eqs (B -58 - Br7l) 1 

i ' ' 1 

pressure/downwash aerodynamic influence coefficient (AIC) matrix 

Laplace transform complex variable 1 

I 

time, sec. 1 

transfer matrix across shock wave, Eq (B-26, D-15); transfer I 
matrices between surface downwash an d sonic line distu rbances,. 

Eq (B-33); transfer matrix between disturbances entering shock 
wave and sonic line disturbances (Eq (B-34)) 

absolute velocity vector 1 

i 1 

local freestream steady flow velocity; also, streamwise velocity 

perturbation (backwash) in shock wave and sonic line transfer . 

functions (Eqs B-12 - B-39) 1 

perturbation velocity normal to shock wave, Eq (B-19) I 

vertical velocity perturbation (downwash) 1 

| 1 

streamwise coordinate with origin; . at wing leading edge 1 


streamwise coordinate with origin . at intersection of sonic 
line and the wing surface 

vertical coordinate with origin at wing chordplane (measured 

perpendicular to wing chord plane) 

1 _____ _____ __ _ ___ 
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-^) relative vertical coordinate for the k^* 1 layer of the 
local supercritical region, with origin at the inter- 
face between the and (k-l) layers. (For k -/l, the 
origin is at the wing surface). See Eqs (A-15 - A-l8, 
C-7) I 

reduced frequency parameter in shock wave transfer matrices, 
Eq (D-2), also coefficient in constraint equation , Eq (E-13) 


Prandtl-Glauert factor, Eq (A-35.) j also constraint 
matrix between freestream and outermost supersonic layer 
Eqs (E-7, E-12) 

see Eqs (E-l6, E-17) 

interface streamline slope, Eq (E-19) 

I ( 12 ) 

difference operator; also AE, ’ ' has a special meaning, 

I 

Eqs (E-24, E-31b); also determinant, Eqs (E-49, E-50) 


C 3 

{} 

FiJ 


dummy variable of integration, Eq (C-13) 
matrices defined by Eqs (B-53 - B-57) 
dummy variable of integration (Eqs A-22, E-85) ) 

3-14159 • • • 

mass density of fluid 
summation operator 
state vectors, Eqs (B-33 j B-4l) 
perturbation velocity potential 

•=2nf ,i"pi rc ular frequency, rad. /sec 

1 

1 

Special Matrix Notation 
rectangular or square matrix. See Appendix B 

column vector matrix (one .column only), See Appendix B 

— — 1 - _ _ _ - 

unit diagonal matrix (off (diagonal elements are zero) 


see Appendix B. 
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SPECIAL PARTIAL DERIVATIVE ROTATION 


= 3 ( )/3 x, etc 
= '3 2 ( )ld z 2 , etc 


OTHER SPECIAL NOTATION 
Laplace transformed quantity, Eq (A-5) 


complex amplitude for case of harmonic time dependence, | 

Eq (A-l - A-4); also in Appendix B, used to define modified^ 

matrices, Eqs (B-53 - B-57, B-60a, B -68 - B-70, B -76 - B- 78 ;) 

i 1 


used to define modified matrix, see Eq B-60b 
vector quantity 

SUPERSCRIPTS AND SUBSCRIPTS 


I 


airfoil or wing surface 
ahead of shock wave | 
aft of .sho.c.k_way_e _ 4 


homogeneous (solution of a differential equation), Eqs (B-l|, 

C-13) 

interference, Eqs (B -58 - B- 65 ); also, imaginary part, Fig.! 
(A-2). 


quantity evaluated at a particular layer interface, Eq (E-5) 


superscript denoting a 


certain interior layer and its assodi- 

I 


ated values of velocity, density, speed of sound and disturbance 
variables, Eqs (A- 19 , C-l, E-5) l 

! 1 

section lift (coefficient), Fig 3 

[ 1 

local (freestream) Eq (l) of text; also lower wing surface,! 

Fig (B-l) , Eqs (B-7, B-13); also refers to layer in multilayer! 

Eq (B-12). ; I 

1 

refers to multilayer, Eq (B-12) i 

section moment coefficient "(Fig 3) ' 


8 






refers to interface between outermost supersonic layer andj 
external freestream I I 


normal component I 

I 

refers to wing surface pressure, Eq (B-42) I 

I 1 

refers t.o particular integral (of a non-homogene ous differ- 
ential equation), Eqs (A-15 - A-17, A-19 - A-22, B-27 - B-06, 
B-42, etc.) J 

real part (of s). Fig A.-2. . j 

1 1 

refers to wing surface, Eq (B-42) j 

refers to backwash velocity, Eq (B-32); also refers to upp,er 

wing surface, or flow field above the wing 1 

I 

refers to downwash velocity, Eq B-32 1 

I 

refers to streamwise direction or backwash velocity I 

I l 

component, Eqs (B-l6 - B-l8) , 


refers to vertical direction or downwash velocity component 
Eqs (B-l6 - B-l8) I 

I 

refers to dependence on perturbation. velocity potential 
(Eqs B-16 - B-l8) 1 

l 

refers to external freestream 1 

I 

I 

NUMERICAL SUPERSCRIPTS AND SUBSCRIPTS . 

I 1 1 

! | 

superscripts refer to |interior and exterior side of a certain 
layer interface, Eqs (lE-1 - E-7)j also subscripts refer to' 
elements of 2 x 2 transfer matrices, Eqs (E-4, E-7))* [ 

I I 

refers to quantities evaluated from purely subcritical flow 
theory, Eqs (B-3 - B-6) I 






MULTIPLE SUBSCRIPTS 'AND SUPERSCRIPTS 
HAVING UNIQUE MEANINGS 


AA,AB,BA,BB partitions of aerodynamic induction and influence coefficient 
matrices Eqs (B-3, B-lo) j 

AB shockwave transfer relationship Eqs (B-26, D-ll) I 


APSH 


p A Tr 
U 


disturbances impinging on upstream side of shock wave caused 

by particular Integra], terms includes (0PSH, UPSH, WPSH), I 

Eqs (B-32, B-34) I 

I 

conditions on downstream side of shock wave i 

I I 

homogeneous, multilayer (Eq B-12) I 


particular integral contribution to multilayer solution, I 
Eq (B-12) | 

refers to wing surface pressure beneath supercritical region 
Eq (B-42) I I 

I I 


shock wave 

sonic line 1 I 

i 

caused by wing surface downwash ahead of shock I 

interference downwash 'induced at lower wing surface ahead Jof 
shock' wave Eq B-45 l 

interference- downwash induced at wing surface aft of the 
shock wave 

I 

downwash contribution | at wing surface ahead of shock wave J j 
due to particular integral term, Eqs (B-32, B-33) I 

| i 

refers to downwash induced by shock wave disturbance at i 

the wing surface ahead of shock wave Eq (B-49) ! 

i ; 

refers to downwash induced by shock wave disturbance at the j 
wing surface aft of the shock wave 


i 




‘BACKGROUND AND PREVIOUS ' RESEARCH 


Review of Existing Unsteady Transonic Flow Methods 
and Exploratory Research 


Older work and its limitations .- It 'is well known (References 8,9 ) thatj 
the equations of unsteady transonic tlow may be linearized for high reduced 

| i 

frequencies; however, it is not clear even in this case that one can neglect 
the sonic line reflections, associated with a supersonic (supercritical) I 

I i 

region embedded in an exterior subsonic flow, and bounded by the sonic line 
and terminal shock (see Figure fl ) . i - I 


OGW : Outgoing 
Waves 


M < 1 -* 

CD 


Sonic Line: x = x (z) 
oL 

^ ICRW (incoming Reflected 

Waves) 

' Shock 

Wave 

M t < 1 
.Li 


Figure . 1 Schematic of .Supercritical Flow 1 

Embedded in a Subsonic External Flow I 

Notwithstanding .neglect of sonic line reflections, Landahl's book (Reference X 
_ . , l I 

9) is devoted almost entirely to the linearized theory of transonic unsteady 

flow, except for Chapter 10, dealing with aileron buzz, which contains some 

j 

.usable results, particularly the "shock-compatibility" relation between [ 
unsteady perturbations on the upstream and downstream side of a shock. 


It is recognized (References 8and 9)' that for thin wings, especially low ; 
aspect ratio wings, linearization at higher reduced frequencies is reason- 

i i 

ably valid. Explicit calculation of non-linear thickness effects for ! j 

/ i i 

unsteady subsonic and supersonic flow past slender bodies by Revell (Refer- 1 

I [ ! 

ence lOand 11) indicate that for slender bodies, non-linear thickness effects 
on dynamic derivatives are rather weak. These results are to be expected^ 
qualitatively because three dimensional relief effects via the continuity 
-equatlonp -require' smaller “excess- flow-field' velocities- to 'circumvent a body] 


FORM 7681-2 




REPORT No. 
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of given lateral width/length^ ratio. Thus, one does not find large local 

variations of Mach numbers. In factf, the first order approximation is the 

| ' 

well known slender wing or slender body theory due to Max Munk and ' 

R. T. Jones (see References 12-14). In this case the lateral lift distribuf 

tion for steady or unsteady flow is. independent of Mach number if I 

/ ^ , 1 

BAR « 1; B =y 1 - M and AR is the aspect ratio of a slender wing. The 

^ * * CO j 

application of slender wing theory to unsteady AICs in the form of Equation 
(5) has been explicitly given by Rodden and Revell (Reference 15)- I 


In connection with unsteady linear theories at M =1, several schemes have 
been devised for computation, using ;the sonic limit for the kernel function 
in the well known procedure of Watkins, Runyan, Woolston, and Cunninghamj; 
(References l6, 17) and also in "transonic box method" procedures (Refers 
ences 9, 15, 18, 19). j ] 

! 1 . 

While the essential importance of the non-linear terms in transonic small 
_dis.t.ur.bance_f.low_theories,_'has. -long. heen_r_e.c.qgniz.e.d .in., st eady flow_circies. _ 
(References 4-9,13,14) little has been done for the unsteady problem. For 
steady flow, Spreiter and Alksne (Reference 5) introduced a notion called 
local linearization in which the non-linear transonic term estimated by I 
successive approximation, yielding equations having the same form as forj 
purely subsonic or supersonic flow, but with the local Mach number varying 
with space. These methods have met with some success in problems of slender 
airfoils and axisymmetric bodies at zero lift near Mach 1. However, the; j 
method is a "simple wave" theory considering only outgoing waves and neglect 
ing incoming wave reflections from the sonic line and from the intermediate 

1 

acoustic impedance changes associated with the gross variation of Mach 

number, density and acoustic speed which occurs between the wing surfacei and 

the undisturbed flow. 1 1 

I I 

I j 

The work of Rubber t and Landahl (Reference 20) also considers non-linear 
steady effects on thick airfoils using the "method of parametric differen- 
tiation" - , ~)Dut'-also fails- to consider- incoming wave reflections-from- the-; — 


sonic— Tinev 
12 





Andrew 'and Stenton (Reference 21) have attempted to apply local lineariza- ' 
tion to unsteady transonic flow, and have done some interesting exploratory 

studies on acoustic ray tracing; however, they confined their studies to! 

I ■ I 

two dimensional propagation in the plane of the wing, whereas, supercritical 
flow effects are most serious fortnick large aspect ratio- wings at moderate 
to high lift coefficients. ! ^ 


In this case the ray tracing(aside from geometric acoustic limitations) would 
be more appropriate in the vertical plane, suggesting a "strip theory" ais a 
start on the general problem. This idea is related to the proposed approach 
which regards the supercritical region as a layered, moving, medium into which 
unsteady disturbances (sound waves) ‘propagate. Acoustic impedance changes 
for such media interfaces ' are given, for example, by Morse, Miles and Riibner 
(References 22 - 24). These impedance changes cause reflections from the 
outer flow to airfoil surface. These reflections occur even in the locally 
subsonic regime; there, however, the impedance changes are quantitatively 

1 i 

weaker , an d cruder approximations ma y suffice. .Befor.e_p.uris.uing _thi s_appr,oa c h 
further, some other recent developments in transonic flow theory must be 1 
discussed. ! t 


Recent research . - Several symposia (References 6, 7) and a recent biblio-j 
graphy (Reference 25) are available describing transonic research since 1950. 
Most of this theoretical and experimental research is directed towards the 
steady flow problems of;, (l) predicting pressures and flow fields past given 
bodies, (2) designing delayed drag rise and "shock free" airfoils. In the J 
case of airfoil design there has been a renewal of interest in the hodograph 
method (References 4, 12, 7 , 26, 27), an indirect, transformation method,' 
using velocity components as independent variables, and yielding a linear 
problem for steady flow, at the expense of boundary condition complications. 
The hodograph method appears to be applicable only to steady, two-dimensional 
flow; therefore, it will not be considered further, except it is noted as a ; 

J 

possible tool for description of the steady flow past a special airfoil • 
shape- whose -unsteady- airloads may be v sought. — — ________ — - - 


l3 “ 




There are two other threads in recent ( transonic research which have led to, 
improved steady flow analytical methods. j 

I l 

o Lax-Wendroff type methods (References 7, 25, 28, 29) | 

I 1 

o Mixed finite difference methods (References 30 - 33) | 


The Lax-Wendroff (forward marching in jtime) (References 7, 25, 28, 29) * 

type of method is actually a transient approach to solving steady flow prop 
blems wherein an initial change (in tne surface boundary condition, for 

| i 

example) is introduced, and the asymptotic (in time) limit of the solution; 
is sought. At each time step, the spatial derivatives are calculated by | 
finite difference methods using data for the previous time step. One knowjs, 
at the previous instant, whether the flow is locally subsonic or supersonic 
and can accordingly adjust the finite difference procedures (e.g., central; 
vs. backwards differences). This method may ultimately be applicable to 1 
oscillatory aerodynamics; however, it is plagued by several problems which; 
would appear to render it impractical 'for flutter application at, the present 

time for the following reasons: I 

I 

(1) Excessive computation time. Reference 28 presents one of the 

most realistic and careful of the Lax-Wendroff schemes for I 

steady flow. It is stated by the authors of Reference 28 to 
“ t I 

require about 2 hours on a large scientific computer for a I 

steady 2D problem. Considering that in flutter analysis a [ 

large number of reduced, frequencies, Mach numbers, angles of 

i I 

attack and wing section shapes, may have to be analyzed, such 

a method would have to toe sparingly applied, perhaps only as; 

a check on simpler methods. 1 

I i 

J j 

(2) The. method has not been demonstrated on oscillatory problems. 

i 

(3) The method shows extreme sensitivity to boundary conditions' 

- - - - • - - - it 

and a tendency towards numerical instability. I 

i 

i 

(4) The physical significance of the "artificial viscosity" (used 

- -to- stab-il-ize-the-calcu-lat-ions-) ; is, as -yet ,-con trover sial-. — j — I 




I I 

With regard to mixed, steady flow finite difference methods (References | 

30-33) these appear faster and more stable for computations than the Lax- 

l 

Wendroff (transient) approach. The results of Steger and Lomax (Reference 
32) and Bailey (Reference 33), based on pioneer work by J. D. Cole and,* 

Krupp and Murman, (Reference 30,31) are impressive. However, no unsteady 

I 1 

applications have yet been developed, though such a scheme might be | 

devised in principle. Again, the prospect for unsteady flow would be 1 
computationally tedious, requiring many flow quantities to be defined, at 
each frequency, and at a large mesh of points, both on and off the body! 

Past subsonic and supersonic flutter methods have all used surface aero r 

, 1 

dynamic singularity distributions, without having to explicitly calculate 
off -body flow field data (excluding component interference problems). 1 
Therefore, any finite difference methods (steady or unsteady) such as | 
References 28 - 33 5 would represent a tremendous escalation of computer 1 
time, even if their extension to oscillatory flow had already been compiet 
ed with demonstrated reliability (which is not likely in the near future). 

1 f- 

The above remarks are not intended ^o discourage development of unsteady, 
finite difference methods, which might well be feasible after a few new 1 

generations of digital computers. On the other hand, it would appear that 

1 I 

in the immediate future, some simpler, useful methods are urgently needed 

I 1 j 

which could shed light on the essential features of transonic unsteady flow 
for flutter analysis purposes, and this is the thrust of the present 
theory. 

! 1 

i 

I ! 

Experimental data. - Ultimately it is desired to compare any theoretical 

j | 

result with experimental data. It is noted that some oscillatory 2D ( 

section data, notably for trailing edge control surfaces, is available in 

| l 

References 37 and 38 , using a pressure measurement system described in ( 

Reference 36 . Reference 35 is noted to contain some oscillatory delta wing 

> 1 I 

pressure data. It is contemplated that the emphasis of the present approach 




than for low aspect ratio thin wings. 
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Therefore, 2D data correlations would be a first objective as a building] 
block in a 2D strip theory as a replacement for' present methods for local 


sections, along the span of a finite 
(AR 4 to 8, say), which are adjacent 
flow. 


wing of moderate to large aspect ratio 
to regions of locally ■ supercritical] 


I 


I 


Physical^Foundation for the Present Approach' . 

I 

Layered medium idealization of the steady supercritical flow past a j 

. I - 

2D section. - We shall outline, in jthe next, section a two r dimensional 

I — - — *• i 

strip theory for 'evaluation of unsteady airload on wing sections at local 
spanwise stations where the steady flow field' has chordwise regions of 

I i> 

locally supercritical flow, such as is shown in Figure -*2 . The steady 
flow field will be approximated as a 


i 


turbances propagate (see Figure '2) 


Vortex 

-Sheets'* 


layered moving medium into which dis- 

i 

I 

I 

I ^ Sonic Line: x = x (z) 

' SL 


M < 1 -> 


U 




y — Shock_Ws 

. M 3> 1 

/ M4 < 1 

/\ M2>M 3 | 

i 

Mj_ >1^ 

\ 



: Figure 2 Layered Medium/Representation of .the Steady 
, _ Transonic Flow Field About, an Airfoil with 
an Enciosed Supercritical Region . 


Empirical substantiation of strip theory in steady flow past large I 

aspect ratio thick wings. - Before 'discussing the 2D strip theory approach 

1 "" " — I | i 

some brief justification will be offered, including an approximate method 


for representing finite span effects 


which is consistent with present state- 


First we 
I 


of-the-art (3D subsonic and supersonic lifting surface theory) 
consider some empirical observations! concerning steady transonic flow past 
wings and airfoils with local supercritical enclosures which help substantij 
_a_te_the .us.efulnes.s _of. .strip _the.ory_ on .finite _wings .which are_ thick’ enough, to 

I ' 

_ha.v_e_s ign if ic.a nt_s.up.e r c r it. ic a l_flow_.effe.cts.. l 
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- During the past several years, a considerable 'amount of experimental and ( 

analytical research on wings and airfoils has been conducted throughout ttie 

world (References 6, 7> 25), along with extensive wing design and develop- 1 

ment studies in the aerospace industry for the wide-body jets. One interest- 

i ; 

ing and useful fact emerges which is described as follows: If one calculates 

_ the surface pressure on an airfoil by I methods which do not account for super- 

/ • , I 

critical regions on the airfoil (any of severa'l will suffice) then, in the 
regions where local Mach numbers are predicted to be subsonic, it is found 
that the theoretical and experimental pressures agree rather well despite] 

I the neglect of supercritical effects. Also it is found that some empirical 

i 

methods (somewhat like Sinnott's method, Reference 40) will fairly adequately 

describe surface pressures in the supercritical region. The present theory 

i 

is an advance, since it also accounts for subsonic interference effects | 
caused by the presence of local supercritical flow regions. 1 

i 

I 

The first point above is illustrated in Figure 3 for an airfoil having local 

1 l 

supercri tic al flow. Sim ilar resul ts ijtave been f ound in wing deve lopment ^ 

^ studies in which subsonic lifting surface theory and 3D non-lifting potential 
l I 

| flow methods have been applied. This indicates that a method based on local 

supercritical flow corrections to a 3D subsonic prediction method has had] 

; ' | 

some empirical foundation in steady flow. Also, many studies have been made 

, I 

to validate simplified wing design procedures based on application of 2D | 

airfoil data. These studies have repeatedly substantiated the nearly correct 

prediction of local chordwise pressures on 3D wings using 2D airfoil data i | 

with appropriate sweep corrections, except in the wing root and tip regions. 


As a final point, it is noted that aerodynamic wing development generally 1 

! i 

leads, by deliberate design, to achieving chordwise pressure distributions 

which are nearly the same along the span, except for practical limitations 

caused by necessary spanwise variations of wing thickness. , 

i 

I 

On the basis of the above considerations, it is suggested that a 2D strip, 


if 






aen? 
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theory for unsteady supercritical flow would be a useful tool for aeroelastic 
analysis. _ I 


Empirical evaluation of 3D flow effects . , 

Prediction of local steady flow environment . - The review of analytical J 
— j methods given in the previous paragraphs indicates that reliable steady super 
critical flow prediction has only recently been achieved for 2D lifting air- 
foils (References 28 - 32) and non-lifting axisymmetric bodies (ReferenceJ 33 j) • 
Even in subsonic steady flows, though methods exist, they often give poor I 
results on practical 3D wing body configurations due to the extreme complex- 
ities of the geometries of wing root fillets, pylons, nacelles, etc., which 
all cause substantial interference. Also viscous and flow separation effects 
are especially difficult to assess for complex geometries. 1 


Despite these 3D analysis difficulties, wing pressure data is generally avail 
able from routine aerodynamic wing deyelopment, which must be accomplished 
-pr-i-or— to -commencement -of— ser-i-ous-fl-utjteii-and-g-ust-ana-l.ys-i-s.;— those-da-ta-a-re- -j- 
also needed to substantiate structural loads before even the strength require 
ments for member sizing can be finalized. Therefore, it is reasonable tol ' 
assume that one has steady flow wing pressure data available to aid in con- 

i I 

structing an unsteady supercritical flow theory in the following ways: 


Locally supercritical regions are mapped, both chordwise and, 

spanwise. 1 

I 

The correct steady flow section lift coefficients and surface 
pressures are available jto "tune" a local 2D description of the 
steady flow above and below the wing. I 


Empirical evaluation of 3D effects in 1 


unsteady flow. - The following approach 


is suggested: 


o It is proposed that unsteady airloads on those wing sections! 
having local supercritical regions be evaluated by the theory 

■_ _ i _ j 

described analytically in 'the next section. Thi's method "includes 
_____ * i f 9— j 


O'JKH E 
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some evaluation of subsonic interference effects aft of the. 

I 

shock wave as discussed in Appendix B. 

| < 

In spanwise regions where the flow is locally subsonic, apply 
the usual subsonic unsteady lifting surface theory for unsteady 
airloads. This is suggested because the methods are approxi 

i 

mately. valid in steady flow as discussed above. 


OUTLINE OF -ANALYTICAL THEORY 


Analytical Assumptions for Layered Medium Analysis of 
Unsteady Flow Past 2D -Sections Having Local Supercritical Regions 


Preliminary remarks .- In the following sections the basic assumptions are 
stated for the analysis of the locally supersonic and locally subsonic 
regions. This will be followed in tjJe section, "Analysis Outline", by an 1 


I outline of the analytical theory, which will be a guide to mathematical 
' [ d'etal'Ts’ c'onta’ine'd - in' the~severad~ Appendices— — ~ — 


Supercritical (locally supersonic) region .- The following is a list of 
highlighted assumptions, key mathematical expressions, and pertinent obser- 
vations : 


The steady flow field is 
methods such as Referenc! 
experimental data. 


4- 


known, by independent calculation, by 
je 32, by semi-empirical methods, or from 


The local supersonic supercritical flow regime can be adequately 

approximated in the( ZydrJrection (perpendicular to the chord) 1 by! 

.. } - t 

a finite number of "z" layers, each having local Mach number, and 

thermodynamic state variables which are, at most, a function ofi 


x, and for a given layer' (Z, 


< Z ^ 


k + 1 


) are a suitable 

' _ . 1 


average of the known continuous distribution in the direction 

I 

perpendicular to the surface. 1 
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The "Z "•''layers are separated by vortex sheets (surfaces of j~ 
tangential velocity discontinuity ) across which occur finite 
changes of steady flow properties' (*an'd, "hence, changes' 1 in-" c 


acoustic impedance) 




The acoustic impedance changes and reflection and transmission 
coefficients across discontinuities are governed by pressure and 
flow direction continuity conditions (Miles, Reference 23) and 
with some modification allow use of previously developed results 
of Reference 23. for reflection and transmission coefficients at 
the discontinuities. Numerical evaluation of the theory of 
Reference 23 has provided layering criteria for the purposes 
of the present study. The methodology of References 23 and 24 
is equivalent and is incorporated in the present theory by a 

I , 1 

simpler approach which employs the Laplace transformation with 

respect to the streamwise variable, x. I 

» 

Each "Z" layer in the supercritical region may, in principle, 
"be - further _ suMivTded~ intb^LoMwise* s ub^irTterva l"s of" "cbhs'tantT 

I . I 

Mach number and fluid state properties separated by discontin- 


uous impedance changes. 


The present theory is further idealized 


by assuming a single average Mach number and thermodynamic state 
' \ i 1 

per \'Z" layer. As justification, it is noted that typical | 

steady state pressure data in the supercritical region for many 

airfoils shows a tendency towards a nearly constant local Machj 

i 

number at a short distance downstream from the sonic line (see 

i - 

Figure 3, and consult References 6, 7? 28 and 32 for examples 

of the near constancy of local Mach number in the supercritical 

region). ' I j 

| 

In any of the supersonic subregions of constant Mach number : | 

(defined by the above mentioned layering assumptions) the 
unsteady velocity potential, and all perturbation quantities 
will be governed by the locally linearized 2D unsteady super- 
sonic theory:., _____ 
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~u = "£l + V0 

J_l 

to - -Pl< U L ^ + |f> 


where, if U^, a^, p^, = U^i/a^ are local fluid properties then:' 


^ <«L * |l 


5* / * 


subject to impedance change boundary conditions at each "z I 

1 i 

layer interface and matching the known airfoil surface normal 
velocity (prescribed downwash) condition. I 

| i 

o Only harmonic motion is considered, therefore any quantity 1 

Q : y is of the form 


Q(t) |= Q (w)e 1 ^ (3) 

o For the~local supersonic' regions,' disturbances on the sonic! - 

t 

line x (z) are regarded as known initial values, which I 
SL | j 

suggests the use of a Laplace transformation on x for solv- 
ing the local supersonic flow field. , 


Subcriti'cal (locally subsonic) region. - 


o The layered medium concept applies in principle to the subs.onic 
region; however, the ra;nge of variation of local flow proper- , 
ties is less than in the supersonic region. Therefore, sub- 
sonic induction effects’ will be calculated from standard sub- J 

j i 

sonic theory as perturbations to the freestream steady flow-, ! 

i ‘ i 

regarded as uniform. > j 

o The principal effects on the subsonic -flow field due to the' 
presence of the embedded supercritical region are described in 
detail in Appendix. B. Briefly they are assumed to be as ! 


f 611'ows : 
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(l) The subsonic doublet strength in the region forward ofj trie 


shock wave is reduced to one-half its usual subsonic 


theoretical value 


I 


because of the loss of the contribute on 


ordinarily made bjj the supercritical region above the ( 
upper wing surface. This doublet strength pertains to 1 


the lift contribution from the subsonic region on the 


pressure side of the wing, opposite the supercritical 
region. 

(2) The subsonic doublet strength is changed, from its usual 

I I 

subsonic value, both forward and aft of the shock wave 

position due to two effects: ' 

| 1 

The above mentioned reduction forward of the shock 

I 

also changes the downwash induced at the airfoil | 
surface aft of the shock wave (compared to usual 1 

subsonic theory). . 


(a) 


P) 


4— 


I l 

Addition al i nterference dow nwash at the airfoil 

surface ^caused by disturbances convected through' 

the shock wave from the supersonic region. These ' 


disturbances 


can be represented by a source distpi- 

1 


bution spread vertically across the shock wave front 

• 1 i 

which terminates the supercritical region. 1 - 

! 11 

(3) Because of the presence of the supercritical region above 

the upper wing forward of the shock, the usual anti- J I 

symmetry of lifting pressure above and below the wing ]is J 

destroyed forward *of the shock. Therefore, the upper' j 
| I | 

surface lift is calculated directly from the local su^er-j 

sonic surface pressure, while the lower surface lift is ! 

I ! 

calculated as one-half of the subsonic doublet loading j 
(including the interference effects described in (l) and j 
(2), above). * . ; 
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(4) It is assumed (and proven in Appendix B) that the modified] 

subsonic doublet lo'ading can be expressed as a linear I 

I i 

combination of airfoil surface downwash due to prescribed 

motion. This provides the requisite "aerodynamic I 

influence coefficients" (lifting pressure increments/unit 

I i 

downwash) in the manner required for aeroelastic analysis 

applications. • 


Outline of Analytical Results 


I 


i 


I 1 , 

Preliminary remarks . - Because of the length of the mathematical derivations 

the reader is referred to the Appendices for all of the mathematical details 

Only key equations will be displayed in the following text to emphasize tlie 

important results. The basic equations and assumptions have been described 

in the preceding paragraphs, or in the Appendices. 1 

I 


-The chosen-form-!- 


Formulation ~of~ the-aerodynamlc— inf Iuence~coeffi'cients~ 

for the aerodynamic influence coefficients (AICs) is of the class (pressure/ 
unit downwash) in matrix algebraic foymat. The matrix formulation of these 

results is contained in Appendix B which is the most important body of ana- 

i I i 

lytical results, since it describes the interactions between the subcritical 
and supercritical regions via shock wjve and sonic line disturbances. The 
reader is urged to study Appendix- B closely to obtain the primary thrust of 
the present analytical formulation, ^ppendix A contains the details of the 
local supersonic solution; however, it can be seen from Appendix B that the 
local supersonic solution is just one of the pieces in the overall analysis, 

. ^ ■ t ' 

and that Appendix B is the "big picture" framework whose details are further! 

I 

delineated by the other Appendices. It is now of interest to highlight the 
AIC matrix formulation of Appendix B.’ 


The desired end result is given by Eqs (B-27 - B74) and (B-79 - B-8l). 

In these equations as shown in Figure Brl,p..53> the superscript, A^, refers 
to the upper wing local supercritical flow regime which is terminated by a 
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shock wave whose locus is x (Z). The. superscript, A T , refers to the I 
1 on * -Li I 

subsonic region on the underside of the wing forward of the chordwise posi- 
tion of the upper surface shock wave. The superscripts, B^, and, B^, refer 
J to the upper and lower wing subsonic regions aft of the upper wing shock I 
wave position. TThe^ downwash values due to surface motion are called 
simply W and W , forward and aft of the shock wave, respectively. The , 

1 .ta Ij ^ I 

present analysis is obviously restricted to a single supercritical region; 

| however, the terms "upper side" and "lower side" could refer to "suction j v 
1 side" and "pressure side", and apply equally well to the case of a singlel 
; •'Supercritical region on the lower side of the wing (at a negative steady J 
I ^tate angle of attack, for example). I I 


From Appendix B, Eqs (B-72- B-74, B-79 - B-8l) the desired (pressure/ 1 
downwash) AICs are given in the following partitioned form, which emphasizes 
the difference in the nature of the aerodynamic coupling between wing down- 
wash control points, depending upon whether they are forward or aft of the 
shock -wave. —The lifting— pressure- for-wa-r-d-of— the .shock— wave— -is— ! 

/ A-\ f f L A U(1 ' - ,,1 

A P T = 1 PI - P ( (M 




The lifting pressure aft of the shock; wave is 


FI ■ l 


j h B u i 

= 4 V - P 




The matrices [%] , are the partitions of the pressure/downwash AICs 
which are required for aeroelastic analysis. Either the lumped aerodynamic , 
forces at structural analysis grid points or the generalized aerodynamic . 
forces, can be obtained from these lifting pressures 'by suitable numerical 
integration schemers '(see 'References - 2, 4, ~10 or 14; for examples) therefore, 
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the pressure AICs provide the necessary data for either, collocation or | 
modal flutter analyses. Because the .downwashes, and W , are allowed j 
arbitrary chordwise distribution, the analysis is applicable to any chord- 


wise deformation . pattern, including chordwise camber changes due to arbi-| 
trary static aeroelastic deformations!. I 


In deriving the above equations several building block matrices are employed. 

The first is the basic subsonic lifting surface solution which is assumedl 

available from any of several "kernel function" or "vortex lattice" type j 

solutions (see References 4, 15-18, 39)* These results are expressed in 'the 

1 1 

following form (see Eqs B-3 to B-6) | I 

TT A -r AA T v ABT If A') I 

W. i| K ' K. if p 1 1 

; W B f= ' K^i ‘k° BB jp° B } m 

A B A 1° 0 J l 0 1 


The subscript o implies purely subsonic flow. In the .subsonic case, the ( 
distinction between the A and B regions is dropped and Eq 8 is written a si 




The desired pressure/downwash AICs ar,e then 


where the inverse Is 




:&] 


A partitioned form is useful in the transonic solution and is written as ' 

J, 


This notation is frequently employed in Appendix B. v ' . 

In the case of purely subcritical flov/, the matrices jMi nd j[%j a _ r i j. j 
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-j identical, by definition. In the case of transonic flow, these matrices I 
differ by an amount which is roughly j proportions], to the chordwise extent 
of the steady supercritical flow region (presumed known) times the difference 

i 1 I 

between upper surface lift for purely supersonic flow and for purely subsonic 

flow. As will be discussed in Appendix A, and below, the local supersonic 

solution is proportional to the Garrick and Rubinow/ solution (References j Ul| 

2, 3, 8, 9) times a layered medium factor. l 

I 
I 

Lift asymmetry . - One of the basic distinctions between the lifting surface 

1 I 

solution for a subsonic flow, vs. a transonic flow with an embedded super- 
critical region, is that the lifting] pressure, forward of the shock is nol 
longer antisymmetric above and belowj the wing for the transonic case. Tihis 
is explained in Appendix B, where a technique is presented for separately 
calculating the upper and lower surface contributions to the lifting pressure. 
The upper surface lifting pressure is derived from the local supersonic I 
solution (discussed, more fully in| Appendices A, C, and E). This includes 
the effe.c.ts_o.f _s.on.ic _line_dis.t.ur.banc.e.s_which_are_regar.ded_a,s init.ial_.c.ondi-^j_ 

tions in the local supersonic solution. These initial conditions are in 1 

| 1 

turn coupled to the subsonic lift distribution by means of subsonic flowi 

field induction matrices which are assumed to be available (at least in I 

principle) from standard subsonic theory. * 


Sonic line disturbances. 


i | 

The sonic line disturbances due to the subsonic' 


wing loading in the regions Ay, By, and By are described in Appendix B in •• 

terms of a discrete set of the unsteady perturbation values for the velocity 

1 ! 

potential, the streamwise velocity (backwash), and the downwash. These ' 

discrete initial values are assumed to be calculated from subsonic theory j 

i ! 

at the sonic line locus x (Z ), where Z. define vertical coordinates of 

olj K. • K. 

i si 

the layer interfaces used in Appendices A, C, and E to describe the local 
supersonic region, A^, as a layered medium. From Appendix B,Eqs (B-13 to 
B-15),the sonic line disturbances are represented as follows for K interface 
layers*' 
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'f U ST.I ■ Kl fr-(p 0 A + « l >L + [<l L ]f p ° B ♦»f} 4’ 

( w S lI ° r d fa * -i p ^ L) L + {a b + w b j <|) 

M [S3 fS(p 0 A i S)}‘ ; p3(p 0 B - .p)J (l) 

r ■ 1 

These matrices have K rows and N_ TT columns, if KL T , is the number of ehordwise 

DW DW 1 

downwash points | 

•' I 

The above equations are summarized more compactly in Eq (B-35) by defining 


a "sonic line state vector" 




{ U SL 


The sonic line state vector defined jby Eq (l6) possesses 3K rows and one 

■. .c.olumn. _ Eqs_(.13_to ,15)-ar.e_the.n summar.ized_in Eq (.B=ii--Sby_using_the_-so.nic_- 

C Xsl^ I 

K . , etc. as partitions 

■ °> A f r SL i, 

1 of what will be called sonic line induction matrices K „ and K „ i! 

defined by 1 * 

.-{“sl} = ( :_ i ( 3 o A + L) } + [S]fo B + 4pB f 

' C | 1 J 

These sonic line induction matrices | K ^ and K ^ are of rank 1 

_ !.Lo,AJ o,B_J | 

3K x N_ TI . It is noted that I p and / p ® "jl are known from the subsonic 
DW 1 o <J X o J . | 


are of rank 


3K x N r 


It is noted that 


solution, Eqs (8 to 12), and have the ranks 1^, x 1 and x 1, where 

I ^ B I 

W TO7 and N ' are respectively the number of downwash control points 
BW A DWg i 

forward, and aft of the shock wave, for the- supercritical case. The J 

"sonic line induction matrices" defined by Eqs (13 to 17) are regarded as 

the second set of known building 'blocks available for manipulation in the 

I i 

theoretical development of Appendix B. I 


OP'HF. F'l 
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The terms designated as 


l and: { ApB } 


are the interference lift disltri 


bution in the subsonic regions in A^ |(below the wing ahead of the shock apd 
aft of the shock, respectively). These terms arise because of the presence 
of supercritical flow in Ay. The factor l/2 denotes the loss of the upper 
wing contribution to doublet strength, compared with the usual subsonic J 
theory, as explained in Appendix B. *The remainder of Appendix B describes 


how these subsonic interference load 
terms of the wing surface downwashes 


increments are derived explicitly in 


W 


] and 


to obtain the desired 


(pressure/downwash) AICs described in the preceding paragraphs. I 

^ i 

Subsonic interference downwash due to shock wave . - In Appendix B, Eqs I 
(B-22 and B-23) there is defined a sHock wave induced interference downwash 
at the airfoil surface, which' is assumed as a reasonable approximation, tb 
satisfy the classical"subsonic kerne] function" type of lifting surface | 

relationship I 


r SH , 
% 


aw: 


SHA 


r: 

, K AA_ 

ap- l - 

1 — 

0 

J — . 

— . — 

l 

h- 

>' 

0 

j 

) take 

into 

account the 



_+ _ 



-(-18). 

I 

(19) 

I 


which relate the kinematic (motional) downwash at the wing surface to the | 
subsonic theoretical lifting pressures Wl and {’.'}■ : 

It is now necessary to relate the interference downwash on the left hand 
side of Eqs (l8 and 19 ) to the disturbances impinging on the shock wave from 
the supersonic side. This relationship is expressed in Appendix B, Eqs ' ! 

(B-20, B-21, and B-46)^. j 


AW 
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Using Eqs (17 to 20) eliminates the sonic line static vector yielding, as 

■ 1 1 


in Eq (1-49) ■ 


A T7 SH > 
AW 

a l 

■ 

M 

1 ATT SH i 

l A V J 

1 


K 


WASH" 


I 


AA 


K, 


WASH ' 


K. 


WBSH 


AB 


WBSH . 

^B 1 J 


•1 ,/ A * L% 
~(P„ + Ap ) 


(p 0 b + ap b ) 



~ WASH" 

g wa 


( \ 



, t 

► + 

1 

_ WBSH 


wi 

A 

J— 

L?WA J 


L’ i 


U21) 


The system of Eq (l6 to 21) is solved simultaneously in Appendix B 
(Eqs B -58 and B- 65 ) to yield the interference loadings in the subsonic 
region in the following form 

I 



AA 


_ L P I 


BA 


AB" 


BB 




B 


I,WA 

P B 
LI,WAJ 


(22 


This gives the interference loading as the sum of two -quantities) (l) : a term 
proportion al to the cla ssical su bson ic doublet loads, via so nic l ine d istur- 
bances, and (2)' a term caused by disturbances generated directly by the J 
wing motion, W , in the locally supersonic region, which then impinge upon 

A I 

the shock wave and are transmitted through the shock wave to cause a source 

distribution on the subsonic side of 1 the shock. The interference loads ' 

I 

defined by Eq (22) are expressed in terms of wing motion downwash by 1 

{ 1 A Bl ! 1 

2 P Q j. Pq r using Eq (10) to obtain the desired AIC ( 

format given by Eqs (5 and T) • \ ! 


One of the key relations in the above sequence is Eq (20), which is Eq 
( B-46b) of Appendix B. The first term multiplying the sonic line state 
vector, is the product of three factors; (l) a subsonic source factor, 

-ir w A,B' 


G.bsh’ ] defining subsonic interference downwash in terms of disturbances 
on the downstream side of the shock wave (see Eq B-45); (2). a shock wave: 

T” baT 

transfer factor T relating shock wave disturbance state vectors on 

~L_bh _J| F - /vpqu“1 

either side of the shock (Eq B-4o) and, (3)- a factor, It , which 

accounts for shock -wave -disturbances - transmitted through the- supercritical 

region from- the- sonic -line. — — — ' — 
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Of the three factors in the first term of Eq(B46a) ,the first subsonic source 

factor has not been written explicitly, but is considered a trivial exercise 

• I 1 

in subsonic theory, and is regarded .as definable, in principle, along with 

all of the other "subsonic induction" matrices. The second factor (the 
shock wave transfer factor) is developed in Appendix D, for normal shock I 
waves; this is merely a convenient statement of results given by Landahl ^ 
(Reference 9> .pp* 110-113)* The third factor, relating to transmission of 
sonic line disturbances through the supercritical region, is a matrix state- 
ment of re suits ’which are developed in Appendices A, .0, and E : , in terms of a 
solution, by Laplace transformation on x, in the layered supercritical region* 
The matrix statement of these resultS| (given in Appendix B in Eqs (B-20, j 

B-21, B-27 - B-3'4, B- 38 , B— 42, - B-46) implies the inverse Laplace transfer- 1 

I i 

mation of -quantities defined in Appendices A, C, and E, by a convolution | 

integral method (Reference 42). Evaluation of the convolution integrals hy 

I i 

numerical integration then lends itself directly to a matrix, formulation. | 

I 

I I 

The seco nd term of Eq .(20) (or Eq B- 4ft) likewise con sists of thre e factory, 
the first two factors being the same ones discussed above. The third factor 
in this case defines the propagation, 
supercritical region, of disturbances 


to the shock wave, through the layered 
generated directly by the upper wing 


will be made concerning the nature of 
will follow an outline of the results 


the required layering; then, there 
presented in Appendices A, C, and Ei 


•Prior -to commencement of the Laplace transform solution presented in this 
report, an assessment was made of the I feasibility of using a multiple 
reflected wave approach, utilizing previous acoustic theory results 
(.References. 22 _to .24.) These, result slcalculate coefficients of -reflection 


surface motion forward of the shock wave. This term is like the Garrick and 
Rubinow solution with an impedance factor.. The matrix statement of this term 

[ : i 

in Appendix B implies the inversion of Laplace transformed quantities derived 
in Appendices A and E. 


Criteria for layering of the local supersonic region-. First some comments 


l 






and transmission for plane acoustic waves, of varying incidence angle, | 
impinging upon a layer interface, across which there is a discontinuous J 
change of the quantities, local stream velocity, local speed of sound and! 

\ I i 

density. Computations were made, using the theory of References 23 and 2,4 

to evaluate the magnitude of the reflection coefficient for single reflec-tioa 

of a plane sound wave, with an arbitrary angle of incidence, (defined as the 

j I 

angle between the wave front and the interface) and striking the interface 

across which a discontinuous change’ in local freestream Mach number is I 

I 

assumed to. occur. It is also assumed that the local density and sound speed 
are related to the local Mach number by isentropic flow equations through- 
out, the supersonic region, which is aj good approximation. I 


The results of the above described calculation show, for a single reflecti'on, 

that the reflection coefficient is equal to about .5 times the change in 1 

I 1 

local Mach number across the layer. |This would suggest that the super- | 

critical region might be approximated by layers of sufficient vertical depth 

1 1 

to o al low a chan g e of local Mac h number of the order .of .2 across ,ea,ch_layer,.._ 

In such a layer model, one would findf reflected waves whose strengths are . 

ten percent of the incident wave strength for each reflection . It can bei 

seen from the analyses in Appendices k, C, and E, that many reflections and 

refractions of waves occur; however, 'they are systematically accounted for 

I ' * 

by the Laplace transform analysis ..which has been employed in these Appendices. 

i 1 

I 

There are two basic reasons for not employing the reflected wave approach! 
directly; (l) the disturbances arisiJg from various points generate cylindri- 
cal waves which can be described as a "bundle of plane waves" of varying ! 

incidence angle; (2) for an arbitrary! steady flow field, dependent on varying 

~ i| 


values of freestream Mach number and angle of attack, the width of the layers 

will vary, as will the sonic line locus and shock wave positions; therefore, | 

1 . : 1 

it would be logically complicated to describe how many wave reflections havej 

contributed to the pressure at any fixed point on the airfoil. 
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;For the first of the stove reasons, the reflection coefficient at an interface, 
due to even a single point source (leading to a Garrick and Rubinow typejof 
disturbance) must be obtained by integrating the reflection coefficient for all 
angles, of - incidence of the plane waj/e bundle which represents the point| 

source (cylindrical wave) disturbance. On the other hand, the layer inter- 

I 

face impedance ratios are shown in Appendices A and E, to be point functions 
of the complex Laplace transform variable which simplifies the handling of 


layer interface boundary conditions. 


The inversion of the Laplace transt 


form of the impedance leads to a convolution integral for the reflection ■, 

I t 

frbm the interface, of a point source disturbance; therefore, this method 
is, in principle, equivalent . to the integration over a bundle of plane | 
waves described above. 

i 
I 

One further subtle point must be made to justify the physical basis for the 

I i 

use of the layered medium model for jthe supercritical region. In the | 

analyses of References 23 and 2b, the impedance ratio across the layer 1 

1 .2 1 

interface is propor tio nal to th e ra tio of p-a^ vaTues where , p, is densit y 

"j ~ 2 - ' 

and, a., is the speed of sound. For a perfect gas, pa- is proportional to 
static pressure, which must be continuous throughout the fluid. However, 
because of the steady flow field streamline curvature, there is a gradient 

to the static pressure in the direction perpendicular to the wing surface. | 

i 1 I 

The present model replaces the continuous static pressure variation by step:-' 


wise changes 


at the layer interfaces which can be regarded as similar to 


membranes capable of supporting static pressure differences due to the steady 

! 1 

flow curvature. However, to the incoming incremental aeroelastic disturj- ■ 
bances, the membrane appears as a porous wall, having continuity of stream-; 
line slope and equality of incremental pressure on either side of the 

; t 

membrane. Thus, the finite layering method chosen here is similar to an: j 

approach sometimes used for example,' in evaluating the effects of atmos- j 

pheric pressure variation upon the propagation of a sonic boom from high; j 
altitude towards the ground. 
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Layered medium- analysis of the local supercritical region . - The following 
is primarily a qualitative discussion) of Appendices A, C, and E. These | 


Appendices contain the detailed mathematical development. 


Appendix A contains a formulation of the boundary value problem for the small 
amplitude, steady or oscillatory disturbances to the local supercritical ( 
region, having a non-uniform steady flow as shown in Figure 1, which is I 
further idealized as a layered medium as shown in Figure 2. The physical 1 
basis of this idealization has- already been discussed. Eqs A-l to A-4 I 
describe, respectively, the convected (wave equation governing the small | 
disturbances, the wing surface downwash boundary condition, and the first I 
layer interface pressure and streamline slope continuity conditions. Eqs [ 

A-5 to A-l4 show these same equations after Laplace transformation on the 1 
streamwise variable x, and taking into account initial values of disturbance 
quantities on the sonic line, which are temporarily regarded as known quanti- 
ties. (Appendix B, in fact, shows how they are coupled to the subsonic J 

loading_to. close_the_analy.tical_feedback_loop.!-)- Eqs -A— 1§— to-A— 18 -show— the' 

analytical solutions within each of two interior supersonic layers, bounded 
by an external subsonic freestream. These solutions indicate the presence 

of particular integrals required to satisfy the non-homogeneous "forcing * 

„ , I ! 

functions appearing on the right hand side of the Laplace transformed version 

! I 

of the convected wave equation (Eqs A-,6, A-7 and A-19). Eqs A-20 to A-23 1 

i i 

show the local supersonic surface pressure as a convolution integral of the 

classical Garrick and Rubinow solution for unsteady disturbances to a stea'dy | 

supersonic stream (References 4l, 2, 3, 8) plus a particular integral contri-j 

bution arising from sonic line disturbances which are induced by the subsonic 

loading. Eqs A-23 to A-26 outline the inversion of the Laplace transforms, 

emphasizing the nature of the integration contours required within the 1 

Laplace transform complex variable (s) plane. Eqs A-27 to A-35 display some 

of the layered medium transfer functions which modify the Garrick and Rubinow r 

solution. Appendix E contains a rigorous derivation of the results given by 1 

' 1 I 

Eqs A-27 to A-35> for the case of two interior supersonic layers bounded by 

an external' freestream': ~ Also contained "in - Appendix - E is 'the" technique “for, I 





J generalization to an. arbitrary number of interior supersonic layers. 

' I 1 

I 

Appendix C contains a rigorous derivation of the particular integrals to j 

_ the non-homogeneous , Laplace transformed, convected wave equations within 

~ the layers (Eq A- 6 , A-7, A-19, C-l). .Besides the general theory, descritjed 

_ by Eq (C-l - C-5, C-13 to C-17), there is also presented an explicit set of 

I 1 

results for the case where the sonic line disturbances vary linearly with, 

Z, the distance perpendicular to the layer. This solution would be a good 

' , 1 

approximation if several layers are used to describe the supercritical | 

I I 

! region, and could influence the choice of layering criteria for a computer 
program development. Also contained^ in Appendix C is a particular integral 
for the external freestream disturbances, wherein, the initial values are 

defined on an extension of the sonic line, which may be arbitrarily located 

I 

for the convenience of the analysis without loss of generality, as long as 
the subsonic induction effects (see Eqs B-13 to B-15) are properly calcul- 
ated. I ! 


-Append'ix - E ‘coh'tTins" a"h e3^ric"iT~derivat - ion of~the~ particular integral contri* 
butions to the airfoil surface pressure which is displayed in Eqs (A-20) land 
Eqs (E-95, E- 96 , and E-97)* The results of Eq E-97 show directly how the 
effects of the layer interface impedance changes influence the sonic line 
disturbance contribution to the airfoil surface pressure in the local super-| 

I 1 i 

sonic region. . , 


EXTENSIONS OF THE PRESENT WORK 


Analytical Effort 


The present theory is completed insofar as the problem is solved, in 
principle, provided the reader accepts -the layered medium model - 
as a. useful mathematical and physical approximation to the total problem: ! 

Appendix B defines the necessary accounting procedure for calculating aero- 
dynamic influence coefficients for aeroelastic analysis applications. It is 
. clear_that_each_matrix_in .Appendix, ^defines _a computer ^subroutine require-.! 
m e n t_,_ and. eac h of t hese subrout ines ma v require subroutines. 1 1 


- . - k e i- 
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It appears, therefore, that the primary analytical requirements will he as 


f ollows : 


1. Write explicit equations for 
defined in Appendix B. 


4 . 


I 


defining each element of each matrix 


I 


Define and write explicit equations for subroutines to the matrices 

of Appendix B as required for the more complicated matrices of J 

Appendix B. I 

I 

Define and write explicit equations' for subroutines to invert the 
Laplace transforms defined by Appendices A, C, D, and E, whose inverse 
Laplace transforms are necessary to define some of the matrices I 
specified in Appendix B. This task might entail several subordinate 
subroutines such as: ! I 


a. A subroutine to find poles of arbitrary transcendental functions 
of the complex variables! This may require a subordinated sub- 
routine to find the zeroes of the complex function defining its 

~ “denominator^ see Eqs - A~27~ to - A“33T* for example )~ “ "j 

b. A subroutine to locate bjanch points in the complex plane for 

an . arbitrary number of interior supersonic layers bounded by a 
subsonic external stream |( see Eqs A-23 to A-26) . I 

| 1 

c. A numerical quadrature subroutine to evaluate the line integrals 

I . i | 

along path segments on either side of branch cuts (line segments 


connecting branch points 


in the complex s plane. See 


I 


Reference 42, Appendix A] Figure A-2 and Eqs A-23 to A-26). 

I i 

d. A subroutine to calculate complex values of residues at the ( 
poles defined by subroutine 3a, above. (See also Eq A-26.) i 

Write an interim analytical summary report covering the above items 
in sufficient detail such that, in conjunction with the present , 
report, a computer programmer could work reasonably independently i 
to begin coding the various subroutines with only occasional clari- , 
fica.tion by_ the. .theoretical .aerodynamic ist, or aeroelastician. • ■ 





Review and monitor computer programming effort; aid in check- 

I 

out of subroutines. I 

Define sample • problems to evaluate the sensitivity of the j 
computer program to the number of layers, integration proce- 
dures, etc. I ^ 

I 

Write a final report according to NASA specifications. . 


Computer Programming Effort j 

I 

Some of the computer programming effort could commence immediately, while 

some of the subroutine coding would . await the detailed definition in the I 

interim report defined as Item 4 under analytical effort. The. programmer 
task activities are visualized as follows: 1 


1. Construct, a flow diagram based on Appendix B of the present 

report, and submit to analyst for review. I 

I I 

2... — Code-and check-out. thej-var-i-ous— subrout-i-ne-s— i-mpli-ed -by ! 

Appendix B and the flow diagram, as soon as sufficient detail 

has been supplied by the analyst in the form of preliminary 

appendices to the interim analytical report. * 

' 1 

3. Upon receipt of the interim analytical report, complete cod- 

! i 

ing and checkout of the various subroutines defined above. ( 

| l 

4. Assemble the entire computer program and write a separate | 

report which would be a user-oriented document referring to 

the analyst's final report for the theory. ( 

| I 

5- Review with NASA the format for the user document, and the 

I 1 

software language compatibility requirements for various 1 
digital computer systems for which NASA may wish to compile 
the software. ^ 4 


6 . Complete the user-oriented report documenting the computer I 

i i 

program. ! 

7 . S ubmit computer program documentation report for__r.e.V-iew_bv_L 

analyst and NASA. 1 3 

..ckhfed' 7 i — • — 3 
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Comments on Matrix Sizes and Computer 
Storage Requirements 


Appendix B defines the matrix accounting system required for solution of the 

’ j I 

] problem. Most . of the . arrays, will be sized either by the number of chordwise 

locations where surface pressure and motion downwash are specified, or byj 

the number of layers required to accurately model the supercritical region. 

It has already been suggested that layers might be defined at vertical posi- 
>•" i ; 

tions having increments of steady flow Mach number of the order of .2 or 

less for greater accuracy. Examination of airfoil data shows that the upper 

i 

j range of local surface Mach number is about 1.4 to 1.6, for an embedded j 
l supercritical region within an exterior subsonic flow. This would suggest 
that 3 to 6 interface layers would certainly suffice, and possibly even the 
two-layer model, derived explicitly in Appendices A and E, might provide ( 
some useful guidance, at least for preliminary calculations. 1 


The number of chordwise downwash control points is governed by two criteria: 

| I 

1. A set of points required by solution of the classical subsonicl 

lifting surface problem. -This could range from as few as two * 
points per chord (for rigik body motion) to 10 or 20 for I 

structural vibration modes having appreciable camber bending'. . 

. It would probably be judicious to define at least 10 chordwisel | 

control points on either slide of the shock wave for the purposes | 

of Eqs (B-3 to B-10). This would define 20 x 20 matrices for a 

strip theory treatment of a wing span station having local super - 

! i 

critical flow. \ 


2. A second set of chordwise bontrol points is needed to develop • ! 

the surface pressure vs downwash relation on the upper wing ! j 
surface beneath the layered supercritical region. This set ofj 
points must be sufficient to allow accurate evaluation of the conp 

i 4 

volution integral described by Eq A-22 and the inverse Laplace j 
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transform of Eq (E-96). The inverse of Eq (.E-96) can be eval-- ; 
uated at an arbitrary number of points by a contour integral • 

similar to Eq (A-23). In evaluating the first term of Eq (A-22),. 
enough chordwise points must be used for the accuracy which is 
required by the numerical integration scheme employed (e.g., 
Gaussian's quadrature, Simpson's rule, etc.). The required values 
of surface downwash at these points can be calculated by inter- | 
polation between the chordwise locations of the downwash control I 

‘ 1 

points used in the subsonic lifting surface solution. This set ! 
of integration control points for the supercritical region might ! 
range from 10 to 50 points depending on the sophistication of the 
quadrature method. In view of the complexity of the logic forj ; 
the computer program as a whole, it may be preferable to begin* 

1 

with a simple quadrature scheme using a larger number of chordwise 
integration, points. I | 

: I I 

This section will.be. concluded by an example., for_ a., refined analytical model,-;., 
of the size of some. of the matrices in Appendix B. It is assumed that five j 

layers will adequately describe the supercritical flow gradients, and that- ! 

l 

it is valid to assume a linear variation with Z of sonic line initial dis- 
turbances across each layer (see Appendix C, Eqs C-6 to C-12). Then as 
indicated in Appendix B, three quantities will be required to define the 
sonic line disturbance state vector for each layer interface (eqs B-13 to 
B-15 and B-35) . For the assumption of ten downwash control points on each' 
side of the shock (a refined description) then Table I defines the rank order; 
of the largest sized matrices appearing in various equations in Appendix B> 
From the pattern in Table I one can easily deduce the ranks of the remaining . 
matrices defined in Appendix B. It appears that the largest individual 
matrix partitions are of the order 10 x 15; therefore, it appears that even 
for a refined layering model, the present method is feasible within the 
storage capacity of many current digital computer systems. This also may 
include even the smaller storage allocations assigned to time-shared, remote 
terminal computer systems, provided the operations sequence -is carefully -- • 
programmed. It seems possible that a respectable" accuracy could'be ^achieved, 
with an active computer core storage capacity of the order of 16 000 words. 
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TABLE I. EXAMPLES OF MATRIX RANK SIZES 

I 

DEFINED BY EQUATIONS IN APPENDIX B 


Assumes 5 layers and 10 downwash 
points on each side of the shock wave 


• C 


Equation 

Number 


B-10,11 

B-12 

B-13 r 15 

~B^TG7T9~ 

B-20r;21 

B-22-:25 

B-26 

B-27-29 

B-30,31 

B-32 

B-33 

B-3U 

B-35-37,39 

B -38 

B-40 


Number 
of Rows 



Number of 
Columns 

/ ■ 
Type 

10 

square 

5 

rectangular/ 

10 

rectangular 

5 

square 

5 

square 

10 

square 

15 

square 

10 

■ square 

1 

column 

15 

rectangular 

15 

rectangular 

15 

square 

1 

column 

10 

square 

15 

square 
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CONCLUDING REMARKS 


A new two-dimensional strip theory for unsteady transonic flow has been I 
presented for the calculation of unsteady aerodynamic influence coefficients 
(AICs). The method is valid for arbitrary, rigid body or aeroelastic defor- 
mations of a wing section whose steady .'flow field is characterized by a lo'cally 
supersonic (supercritical) flow region which is adjacent to one side of thje 
wing, and which is embedded in an exterior subsonic flow field. In the I 


absence of supercritical flow effects (the present theory reduces automatically 

to classical subsonic theory, which is desirable. I 

. I 

I 1 

Based on empirical evidence for steady transonic supercritical flow past air- 
foils and wings, it is believed that the present theory, when fully developed 
for routine digital computation;, will provide a valuable tool for the aero'- 

i 

elastician. It will be suitable for performing aeroelastic calculations for 
-moderately- swept.,— large aspect-ratio, -j-thick-wings -oper.ating-at -high-sub sonic- - 
speeds, and will use simple sweep theory concepts to define the necessary j 
equivalent two-dimensional freestream flow properties along the swept chord 
direction. * 


The present report also tentatively suggests empirical means for estimatin'g 

| i 

finite span effects, using existing subsonic lifting surface theory, plus ( 

experimental^ data for wing pressures in steady flow. The suggested techni 1 - 

| I 

que defines the chordwise and spanwise extent of supercritical flow from the 
measured steady flow wing pressures, and replaces the subsonic air loads by 
the present new theory for those spanwise stations sharing chordwise localjly 
supercritical flow. 1 


(.Various possible theoretical approaches in the literature were reviewed and J 
! rejected on various grounds, prior to ^selecting the present layered medium j 
j theory, which is believed to provide a middle ground in computational complex- 

! ( i ! 

• it y, “and" ha's the "advantage "of "providing "go'od" physical" insight - into"the~ — 
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T 


crucial aspects of the problem. Merely by formulating the analysis by the 
present method yields (Jmpbrtant iconclusions^rs * For example, the following 
conclusions have been drawn from the present studies: j 

\ U 1 , 

1. The usual anti-symmetry of lift above and below the wing is altered 

to an asymmetry of the upper and lower surface lift forward of the 
shock wave. These contributions must be separately calculated. J 

2. The air loads adjacent to locally subsonic regions can be calcuj- 
lated by the usual doublet procedure, if the air loads ahead of I the 
shocks on the locally subsonic side of the wing. are. .assigned . acfactor 
of one-half, and if additionally one accounts for subsonic downwasli 
interference: at' the wing surface caused by disturbances transmitted 


through the shock wave from 


the supercritical region. 


3. If one defines suitable subsonic flow field induction matrices,! 
based on existing theory, one can couple the subcritical and . 
supercritical regions via shock wave and sonic line disturbances 

— wKiclT - are then capable . of algebraic elimination by their expression 

I . 1 ' 

as linear combinations of wing surface downwash; hence, it is I 

possible to explicitly derij/e the desired AICs required for aero- 

elastic application. . 

l 1 

4. A preliminary assessment ofj the computational aspects indicates] 

only a modest computer storage requirement for the present theory. 

| l 

The use of two to five interior supersonic layers may well suffice 

for many transonic flow problems of interest, where local wing I 

surface Mach numbers seldom) exceed 1.6. . 


5. The layering effects in the 


supercritical region are significant, 


as the impedance change between the wing surface and the freestream 
is proportional to the ratijj of local to freestream absolute static 
pressure. This can vary- by 1 a factor of 3 Tor a typical local Mach 
number variation in the flow field between 0.75 and 1.6. 

| " ' ] 

6. Unpublished numerical studies by the author, using acoustic methods 

have shown that for a wide range of conditions, a Mach number change 
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single reflection from an interface in the local supersonic region. 
This result provides a simple ^physical basis for defining the j 
layering criteria for a given steady supercritical flow field. I 


7 . The nature of the local unsteady supersonic solution is easily 1 

! 

understood in the present'. theory, wherein it is expressed in the Iform 
of the classical Garrick and Rubinow solution for unsteady flow dis- 
turbances to a uniform supersonic stream, times a layered medium , 1 

multiple reflection factor, plus a particular integral term arising 

I ' 

from sonic line disturbances. The sonic line disturbances are induc- 
ed by ' the: wing surface air 'loads adjacent to the subsonic region. 

» I 


I I 

'It is strongly urged that the present theory be implemented for digital compu- 
tation, since it- is inherently well oriented towards aeroelastic analysis' 1 
reauirements bv the naturp of thp formiil ati on . ' 
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APPENDIX A 


LAPLACE TRANSFORM SOLUTION TO THE FLOW FIELD IN 
THE MULTI-LAYERED, LOCALLY SUPERSONIC REGION A^ 


The local supersonic region may he broken into layers of constant local , 
Mach number with' increments of about .2 in AM according to the results 
of studies alr eady mentioned. The corivected wave equation in each 

region is solved by Laplace transform on x, assuming initial values on tile 
sonic line are given. A solution will be described for two interior super- 
sonic layers I Layer 1.'. is defined by(© ^ Z ^-Z^) and layer 2 is defined by 
(z 1 ^z<z 2 ). I 

The velocity potential equation for harmonic '.motion is (omitting e 1 ^): j 


0. = B* 0-, + 2i«A M, 0 1Y - JSL 

-Izz- — 1 hx ~ ' a; aqf— 4 - 


- (A=ll 


The surface boundary condition is 


0 , (X,0) =W (X); W (X) =W (X,t)e“ 

h a - a a 


( A- 2 )i 


The pressure continuity at Z = Z^ is: 


P-L (X,Z 1 ) = p 2 (X, Z 1 ) 


- p i (u x 0 lx + m'-p = -p 2 (u 2 0 2X + 


The streamline slope continuity at Z = Z^ is: 


§• h < x > V = I *2 (x - Z U 

U 1 X Z i 2 Z 


(A- 3 ); 


(A-4); 


' ~KHfrE O 
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Where 0^ and 0^ are the velocity perturbation potentials in layers 1 and 1 


separated by an interface at Z = Z, (see Fig. A-i) 


Extended! 
Sonic Line 


Sonic | 
Line J 


'J u »1 

/ ^~ ' \ Z 2 


M = U /a <1 


CO 05 00 


Shock 

.Wave 


“2 = ^V 3 ^ 14 ^ 1 ) 
M 1 = ^/V 3 !^! 2-1 ) 


A 1 I 1 

Fig. A-l Two Supersonic' Layer Model Bounded by an , 

External Subsonic Free Stream l 

i 1 

I 

We next Laplace transform Equations [(A1-A4):' with .respect to x and obtain 


? <S,Z) . 0 (x,z) ax 


(A-5) i 


ft - K A = 3 1 (S ’ Z) 

ZZ | 


{■(A-6) 


ft -K ft = G (s,Z) 
ZZ 1 


(A-7) 


where 


G x (s,Z) 


B I 2 < s *SL 1> « 


— “1 SL 
a l 


G 2 (s,Z) = - }b 2 (S0 s ^ 2) + 0 X 8 [ Z) > + 2iU L M 2 0 s f } [ 

V •- • a 2 . 


= b/s 2 + [j^i^s - iilj • 

I L 1 a, J 


(A-8) 


(A-9) 


(A-10)| 


Kg = Bg s + 


[2M 2 i(-f-)s - £3 J 


2 a 2 ■ 


(A- 11) 


Here •0g 3 ^-,-0 x ^g'^, _ 0g 2 ^ , ‘0£^ 2 ^ ~are “the sonic line ('initial) values o'f 
velocity-potential— and - streamwrse _ perturbati'on _ velocity a.t the sonic 
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line for layers 1 and 2, respectively. 


The boundary conditions are Laplace transformed as follows: At the airfoil 


surface 


X (s,0) = W (s) 


X Z ' a 

At the first interface, pressure continuity yields 


(A-12), 


P x (UjS + iO)) % (s,Z^)= P 2 (U 2 s -!- ill)) 0^ (s,Z^) 

At the first interface, streamline slope equality yields 


k\ (s>Zi) <s>Zi> 


. (A-14). 


Solutions for $ (s,Z) and '$ 0 (s,Z) satisfying Equations(A6) and (A7)can. 


he written 


~$j~ — (-S-J-Z-) = -cosh— (-KjZ ) — + -Dj—s-i-nh— -zl) (A-15) 


z = Z for (0 < Z < Z ) 


^ 2 (s,z) = C 2 cosh (K 2 j z) + D 2 sinh ( K^z) + (s,z) 


(A-16) 


z = Z 


- Z ] _ for (Z < Z < Z 2 ) 


Solutions of this type are written for each layer, except the disturbance 
to the free stream is of the form ' I 


(s,z) = C ro e _K e° Z7 + 0 p . (s,z) 


(A-17) 


2 = 2 z 2 


(A-18) 
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If is subsonic, then B»=+ i /i-M and obeys rules of complex 
variable theory. We imagine the initial conditions given (by subsonic 


line (see Fig. A-l.). The 1 
, are solved by applying the | 
and(Al4) at the interfaces and by I 


theory) on an extension of the sonic 
coefficients C^, D JL D 2’ , * ,C N ’ D N j 
boundary conditions of the form(Al3) 
satisfying the surface downwash conditions. 

The solutions, {$_ |, are particular integrals satisfying the non- 

^ F , | 

homogeneous RHS of Equations ( a 6,A7> etc.). 


C^,h) 77 - ?k 2 *p,k (sjZ) = \ (s ’ z) 

Z1Z1 * 


(A-19) 


G, (s,Z) can probably be regarded as a constant or linear in Z over a 


narrow layer (Z, . < Z < Z^) based on] 
the sonic line from subsonic theory; 
easily be found. 


a curve fit to induced velocities on 

I 

therefore, solutions to(A19) can 1 


I 

4 - 


The Laplace transform for the upper -surface lifting pressure can be 
written in the form 


Ap- (s ,0 ) = E 21 (s) ftp; (s,0) + AP Fjl (s,0-) 


GR/ 


where (Ap (s,0) is the Garrick and Rubinow solution 

a ^GR 


Ap (s ,0 ) = -P 1 (U^s + ixu) R 


Wa(s) 

1 


1 

(A-20) ' 

♦ - 

I 


(A-21) , 


and ? A|?p (s,0 ) is the additional surface pressure caused by sonic line 1 


disturbances , 


The oscillatory surface pressure is obtained by the inverse Laplace 
transform of (A20),. i 


' I 

1 

1 i 


GR 


A? (x£>) - E 21 (x-§) Ap (?,0) d? + Ap p 1 (x,0 ) 


(A-22) 


- T 

d* 



Using the convolution principle (Ref ( . 42) , the function is given 

by the Laplace transform-inversion— theorem ' * 


I /■ 0 + 1“ 


E 01 (x-§) =— I i,E pl (s)e 




c - i“ 


The inversion .integral (A-23)is evaluated, along the indented contour 


shown in Fig. A -2 
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FIGURE A-2 Complex Plane and Inversion Integral' 
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Branch Points: 


S k,l’ S k,2^ 


Defined By': i /, \ = 0 = 


It can be seen from Equations (A10 &" '.All) that K, and K 0 are quadratics 

1 2 l 

in s; therefore, one can write I , 


K 1 = l 7 (S - S 1,1> (S - S 1,2> 


*2 = ± 


■ k h c e n 





I 


and similar expressions for each layer. It turns out that the branch | 
points (s = s^. li- e on imaginary axis for all layers which are 
locally supersonic, and for locally subsonic layers, the branch points | 
are symmetric about .the! imaginary axis . I 


The integral(A23) may be evaluated by the calculus of residues (see Ref.j+2) 
to yield 


E 


21 


(x-§) = 


2rri 


r 

2rri 


N 

E 

J=1 


e sj (x 5 ^Residues ^E- (s.) 

21 - j 


(s)ds fl 



(A-2 6 ) 


The last sum i'n(A26) represents the contribution of all the line integrals 


of if ^(s) over various path segments! shown in Fig.Ar2. It can be shown 1 
that Egj. (s ) is sin gle valued ev ery where in th e com plex 1 s" plane, except on 

li 


opposite sides of the branch cuts (1'ine segments connecting branch points 


sV a ;ahd' .s kP for each layer), 


*-?:k2 


It also turns out that the integrals 


aroundTsmall circles surrounding the branch points vanish, and E 2 ^(s) is| 


an even function of (or in any’ layer separating adjacent supersonic 
regions k and k-l). Therefore, the 


only line integral segment contribu-| 
jfirst layer and the last layer’ 

This has been rigorously proven i 


tions to Equation(A2.6) are from the 
adjacent to the subsonic free stream 

for the special case of two supersonic layers bounded by an. external subj 
sonic free stream. I 


The first sum in Equation(A26) represents.- residues at poles of (s) i 
with contributions from each layer. E , (s) can be written in the form 1 
( numerator /denominator) and the denominator has -the form j 


D e (s) = { 1 - H 21 (s) tanh (K^)} 

- - 1 .4 


(A-27). 
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The required poles are zeroes 'of Equat;ion(A27) which are countably infinit’e. 

~ / N 1 

Contributions from outer layers are contained in 'H (s). I 


It is thus seen that the Laplace transform solution is equivalent to a I 

multiple reflection solution but provides a systematic computation procedure. 

I 1 

I 

Special Solution for Two .Supersoliic Layers Bounded 

by an External iSubsonic Freestream ! 


In this case 


E 21 < b > = N E 


(s) 


(21 


D I (s) 
E | 2.1 


where (s) is given by Equation(A27) and 

^21 


‘ H ir< T ) _ = z 2i 


N. 


21 


(S) = {‘H; 


'("s')~E~ ( s ) 


21 


(s) -■ tanh 


(K,Z,)} 


Do ,(U 2 s + iw) U„K, 

Z* (s) = ^ 

P y |(U 1S + id) ) U^K ^ 


E 2 (s) = 


{z r COj 2 + tanh (Kg (Z^-Z^)} 

\l +%, t2 tanh ( k 2 ‘ ( Z 2 " Z 1 ^} 


and 


't ' D„ (Uo=S> + iU() UeoKp 

, 2 


0 2 (UgS 

The quantities zf 2 i(s ) and 2&,,2.( s ) are, 


+ 


id) UpKco 


(A-28), 
I 
i 

I 
I 
I 

(A-29) 

(A-3Q) 

I 

(A-31) 

(A-32) 

1 j 

I { 

(A-33) 


in fact, the Laplace transforms of i 
the impedance ratios across the first and second interfaces. In evaluation 
of Koa for the case, M» < 1, one substitutes into Equations like(AlO) and | 
(Ali)and obtains: j - 1 
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APPENDIX B 

OVERALL SUBSONIC/TRANSONIC LIFTING SURFACE SOLUTIONS 

I 

Let the subscripts A and (B„, B t , A_ ) denote, respectively, quantities 

U U L L| 

associated with the local supersonic and subsonic regions, respectively 
(see Figure B-l). 



'Figure B-l Designation of Regions of Locally 
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Where H denotes the homogeneous soluti'on (due to the prescribed surface 1 

downwash) and P denotes the part.ic.ular) integral (due to sonic line I 

disturbances, as discussed in Appendix A) . 0 Au and p A U include ._ j 

all of the multi-layer reflection effejcts discussed previously. I 

I 

We next describe the surface downwash 'in the chordwise regions adjacent | 
to A and B by column vector matrices ’ and fe}- It is assumed that one 
has available a- suitable subsonic kernel function/vortex lattice collocation] 
type solution for the lifting surface pressure of the form I 

AA j 


1 % 

W B 


K 


K 


BA 


K 


K 


AB 


BB 


B 


(B-3) 


In the purely subsonic case, the A and) 


B distinction is omitted: 


The- s o-l-ut ion -f or 


; where 


r-{p.o} .i. n -- 


r l r 1 1 

r i- 

i w ; = L , K o J 

i.Po J 


(B-4) 


r i 

terms-o-f-'|W-j 4x-s-:- 


I 

1 J[ 




Q, 


■p J 
^o i 


u K o 




-1 


(B-5) 


(B-6) 


The object of the ensuing analysis is jto replace the upper surface solution J 

’' : in region A,, by the results of the local supersonic solution and modify the 
v | l 

local subsonic lift distribution to acpount for the interference in Regio^ B 

| and A t caused by the presence of the local supersonic flow in Region A... ! 

h'-. I U i 


As a first step, define the lifting pressure as the difference between upper 

I i 

and lower surface pressure: j , 


P = (P 1 L - P D ) 


(B-.7) 


l- J 
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In purely subsonic flow, 




r~Pr 


and 


2p. 


h> 


" 2p tL 


I 

(B-8) 

I 

I 

(B-9) 


When supercritical flow is present in Region A^, this anti -symmetry property 

of the lifting pressure is destroyed, and one must calculate separately the 

| - I 

upper and lower surface lift contributions ahead of the shock. Let us call 

and the lower wing regions forward and aft of the chordwise location* 

of the upper wing surface shockwave location (see Fig. B-1-) . i 


It is still possible to calculate the lifting pressure in the subsonic , 
regions B T _and B by assuming an anti-symmetric lifting doublet distribution; 

; iJJ ■ L I 

i however , "the doublet strength will be changed from its original subsonic ( 

-j theoretical va-l-ues— (-required— ' by-def i-nisb-ion— 1 bo-match the-mot-i-on downwash-)-^-4 — 


Three distinct sources of interference 


downwash arise: 


(1) A change in doublet strength' on the subsonic side in region A^. | 

(2) . The loss of one-half of the original subsonic doublet strength * 

contributed by the upper side of the wing in the (now supersonic,) 


region called A^. 


(3) 


The addition of an interference downwash along the airfoil surface 
in A^, B l , and B y due to disturbances transmitted through the ' 
shock wave. These disturbances can be regarded as caused by a , 


monopole source distribution 1 2 distributed vertically along the 


1 


shock wave front whose strength is everywhere equal to the com- , 
ponent of velocity normal to^ the shock front on the downstream ! 
side (allowing for the possibility of a curved or oblique shock J 
wave ) . I 

The strength of this velocity is given by Landahl's (Ref. 4) shoick j 
-wave compatibility relations., .suitably .expressed in .terms, of _t he'' 
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I I 

local supercritical solution (given by Eq. (EL)) and evaluated on] 
the locus of the shock, (X (Z)). on the upstream side. J 


The final form of the lifting pressure solution will be similar to Eq. (B5) 


in format : 


{ } - M ( w a) + .My 

M . [ s «] {»*} ♦ MW 


(B-10) 


(B-ii) 


t 1 1 

it is necessary now to describe the matrices' in more 1 fundamental terms. From 

Eq.(B2) one can write- the supersonic region surface pressure as: i 

I 

y '[c]ki*ri ( u S l) + mKl] <B ;- i2) 

where ( U-A , / W„A are the X and Z velocities on the sonic line , 

V SLJ l SLJ f r- pM Q ,r- pML -T 

( at 'each -layer -centerline-)-due-to s ubjs onic-dis turb anc e s-and— — tj 

are the induction matrices for the effect on surface pressure due to the | 
particular integrals arising from thej local multi-layer supersonic solution. 
The term ^SLjl, recognizes phase differences in velocity potential 

along the sonic line. 1 I 

i 

The sonic line disturbance potential jand velocities are expressed in terms i 


of subsonic induction effects 


{ 4 - [-c] ( * ^ + + k 


H - [<f] { * e * * ki L J k ♦ y 

( 4 - MJ { m k + k 


(B-13)| 

I i 


(B-14) 


(B-15) 


In Eqs.(B13) and(Bl4) -^Ap Ap^j represent the interference loading in 
the subsonic local flow regions A^, B^, B p due to the presence of locally, 
supersonic flow in region [ AJ on the upper side of the wing. The factor 
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l/2 multiplied by| p Q ^ + Ap^} accounts for the local supersonic region | 
,Ay f not communicating with the subsonic region (except via the shock waye 

s"ource distribution). j 

% * 

I 

Next consider the disturbances transmitted from the supersonic region A- 1 

U I 

through the shock wave I 


<*f>j 

{ Of- [«’“] { Cl ♦ - I<’ S M CB-i 

fe* SH ] {-sh A u } - fe^ SH ]KH u {-fe’ s S^s^l <b.; 

Equations (a.6)’ to, (Bl8) represent a generalized form of Landahl's shock 1 
compatibility relations for curved add oblique shocks. For a given shock 
geometry the normal velocity on the d'ownstream side is I 


(B-I 7 O 


(B-18) 


] . T %)_ F SH f TT %' , £ SH \ L % 

| V SH,n J I n x J { U SH j ! n zJ { W SH _ 


(B-'19) 


where- . n c 
front . 


k n is the unit normal vector to the shock 
z I 


The surface downwash in the subsonic 'regions induced by the normal velocity 
(source distribution) on the downstrejam side of the shock can be expressed 
in terms of suitable Green's functions: 1 

d ] ■ [1 {•,:"] • [wdl - did] '-- 11 




p 
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| V’} =I g s“aJ{ w a + <} + ] U P>] 


w. 




SH 


r ww i c 

L G sh,aJ { 


W A + : AW a ( + HW. 


) TiJ A U V 
P,SH Jl 


where fromEqs.(Bl) and(B2) one can define 


I u < 4 < z » 


f 

J 

<! 

* 

U SH 

= t 

. fc 

J 

-L 

■f 


I 

-V 

„ 

W i 
SH J 

f = l 


30 


*u 


92 ( X SH^'> 


Next, hy analogy -to Eq.(B12) one can 
as a linear coml 
matrix relation 


(B-28) 

l 

I 

I 

(b-29) 

I 


(B-30) 


(B-31) 

I 


as a linear combination of { tjL’ U SL' 


I ) ATT P M AlJ A U | 

^express J A ^ A » 0 p,SH’ W P,SH.’; ,U p,SH ; / 
g L J. . Let us define a partitioned 


Y > 

P 

AW a 

\ 

A u 

0 P,SH 

1 


u * 

P,SH 

> = 

Atj 

w u 
L p 3 sh J 



WAP 

0SL 

c « 

0PSH 

0SL 


WAP 
pSL 

'0PSH 

l USL 


UPSH UPSH 


0SL 


WPSH 

-0SL 


jUSL 


T-, 


WPSH 


iUSL 


WAP 

WSL 

0PSH 

WSL 

UPSH 

l WSL 

WPSH 

L WSL 


t 

— ’V~ 

( 0 

; SL ; 

1 

1 

'< 

; U SL > 


. W BL. 


(B- 


32) 


It is convenient to introduce shock front and sonic line "state vectors" i 


such that: 


Kl -MM 


(b-33) 


(B-34) 
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where the "state vectors " are defined as 




/ W SL l 


; si ?,sh , 

If u Au U 

fp,SH =| P ’ SH f 
P,SH ' 


Eqs.(B27) to(B29) can be regarded as |defining a total shock entry "state 
vector" 

f A, A 

r Aui \ / 

k H -=fsH u < B - 

> J I . At 


From Eqs.(B27) through(B37) 


r A u/I T v H l f 1 T wap! J 1 T 

4°BHy- G WA + 1 T Sl! ’ °SL ’ + ^ ! 

A J l— Jw / -L -J a J - L 


t CT aT 

SL SL 


One can rewrite Eq.(B26) in terms of a downstream shock "state vector" 


t-1 ■ 


■ W 1 


Eq.(B26) then implies 


M 


i _ 

f T BA 1 i L ^ J 

|| SH J j J SH Jr 
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Equations (B13) , (Bl4) , and(B15) can be] • rewritten in terms of a relationship 


between the sonic line "state vector"! 


and the subsonic doublet loading 


■ { ° S l) = [ K o,f| + ^ AL >} + [ K o j L ]fo + 4 P B ] 4) 

- | '-I J \ 

One can by analogy to Eqs . (B2) , (B27) 'through(B£9) calculate the upper surface 
pressure in region Ay in terms of! the sonic line "state vector" 


A tt /I f Att/1 ,I = P - 


p Vo) 




WA 


1 W A + %' 

r .1 


') | Ay| 

/ + p P P j s f 


I 

(B-42) 
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By analogy to Eq..(S32), the surface pressure due to the particular integral 
in 'the supersonic region may be written: [ 

i I 


Au 


\ r 


■P,S 


PV P 


T; 


SL 




Eqs.(B42) and(B43) combine to yield: | 

I r — -v— -r=i- 


U 


P (Z=0) 


Ps^Hj 

°WA 


'W 


* 


A 


T 


,WAP 
SL . 


i 


, f) a SLt ,+ 


pAu j P 

C SL 


Y- 


SL 


Equations(B20) and(B2l) can be rewritten as 


(B-^3) 


(b-44) 


SH 'J 


WAp 


*1? 

'“v 


g bsh 


' SH 1 

► — 

wb! 

5 

°SH 

aw b | 

\ -L» J 


G ^ 

L BSH 

! 


i. 


■l • 


From Eqs. (B45) j(B40) , and(B38) 


1,4 


* - 


B-, 


G. 


wa l1 


BSH 

wb t ' 

L 

J BSH 


c ba 

SH 


r< ^][ T “']l' [ T 


APSH 


(b-45) 


1 - •*. 


SL 


) r ¥ H 1f 

t sl| + [Sja w . 


r (to 


r wa 


■L ■ 


SL 

"'WBy" 

_lt g sl- 


'4~°SL { 


t— 


WASH 
^A 
"WM 

b 5 ^ d 




. 1 
l A j 


(B-46b) 
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in terms 


Substitution of Eq. (B4l) ^ Og^j- in Eq.(B44) yields '^p ^(Z=0)^< 
of ^ ?(p o A + ^p AL ) Jj 5 ^Ap B J ’ Substitution of Eq.(B4l) for a gL || 
J (B-46b) yields ^ AW A ^ H , AW B ^ J as a linear function of £^(p o A + Ap AL )|. 

Rl 


in 

i 


Thus 


AW a ! 

AW. 


SH 

iL_ 

SH 




B t 


WAt 


G 


w: 

"SL 


SL 

o;a 


'1 / A , a a l » 

:|(P 0 + AP ) 


K 


SL| 
o ,B 


B 


A 

+ Ap v 


IV 

— 1 _ * 

= WASH •-] 

Ra 

•R 

t 

, ’ WBSH 

1 {Vf 

/ 

lsja J 



(B-47) 


Substituting Eq.(B47) into Eqs.(B24) and(B25) leads to a well posed set of 

linear matrix algebraic equations for l{ i(P 0 A + ^ and -^Ap Pj}J expressed 

as functions of / |p A } and {p f l . These results may then be back ' 

A r Ay "\ i ° /' | 

substituted to get < p (Z=0)K the upper surface supersonic pressure i 
as a linear function of / W A , |p Q ' , P Q J*. Eq. (b 47) may be expressed 
more compactly as * . I 


X T SH 1 

K 

» — 

„ WASH 

JR A 

AA 

-1 

\| 

^ WASH T 

KJ 

1. 

v WBSH 

I ^A 
1 

i 

T-r WBSH 

Sb 


1/ A ■ A A L'j 
?(P 0 + A P ! 


Ap B | 


^c B + 


J I— 


„WASH 

y a _ 

WBSH 

^A 
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where 


-J 


rc*] 

i 

■ im^i 

i 
i 


r ^ h ] 

LC 3 ”] 


(B-50) 
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Putting Eq.(B49) into(B24) and(B25) yields 

I 


r\Z WAsilJw A A A lF FwASrfl f B . b! F aa] fx &ashI (i 1 
LL K AA J^o + 4 P )t + _,L K ABli P o + - Ap l tL K o J Fo ]fwA J1 Ta| 

; A ] { • } K o AB | { fipB |" ( B-151 ) 

jM{*( p o A * * pAL >} +br% o B * >4 

* MM <*-n 

I I 

The quantities ( P o A } and K B } are! known from the subsonic solution J 
(Eqs.(B3) to(B6)); therefore, Eqs . (B51') and(B52) may be written in the 1 


• < 1 

(B-52) 


.tion 

the 


following form: 


^M + HM = M M 1 Fj MlF (B ‘' 53> 

HffTWM - "H i 1 f pfl IBFF} ^ 


(B- 53 ) 


where 


jRi - M " ■ iDM 


[Fab] = i 

Fi - ; 

Fr H ]J 

1 1 

1 

1 


[vfj* j 

tri 

(B- 55) 1 

IM - 1 
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1 

1 

Ia '“3 - 
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■ 1 ! 
1 I 

i 

Fba] - i 

F BA 1 - j 

Frl 

1 

i 

t | 

, 1 

I he]! = J 

a b 4 - J 

FrJ 

(B- 56 ) 
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Eqs.(B53) and(B5^) may "be more compacJly written as: 


ja_p a 5 

Ap B \ 


1 7 1 at .fcL WASH " 

A'. JJpoil + 5™ _ J 

n i udou 1 


-L J i “ p ~ :l _l ji p o B i [Lv:j i 1 

Solving Eq.(B57) yields finally the subsonic local flow interference 
pressure increment due to the presence of the supercritical flow region 
Ay on the upper side of the wing. I 


iWBSH |! W A 


i Ap 


,WASH 


A +; Ji* V == -- ■ C _°_ 

B| 1 WBSH ! p B B 

. . p o / . g wa _ X h L L P i _ p o .j 


I i A 


° l + _I_jWA 4 W 

”B P B 
o j I,WAj 


'The interference pressure is s"een to b ( e proportional to the basic subsonic* 
-loading_computed-by_class.ical-subsonicj_kernel_f.unction~or_vor.tex. lattice — ^ 
theory, -pluS'-'an' increment from the local supersonic solution. ( 

Local Supersonic. Upper Surface ‘Pres sure . i 


From Eqs.(B4.4),(B4l) and(B5&) one can write the local supersonic pressure ( 
in the following form: 1 

fr] I 


SL~| f 

1 k °> b Ji 


I B . . B 

Po + Ap I 


Letting 


' T WAP. 
. SL 


+ T r 


Ps A U p 


? sAu ] = M iPt wA pi| + f T PsAuP l 

_ G WA J _ L?WA J jl_ T SL J |_ T SL 

z fs\~\ _ T p sVT + UM r sl] r a i 

_^W A J L WA £WA J jiSojAJ [fl,WAj 
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Lifting Pressure Aft of the Shock Wave 


I O I 

From Eqs.(B3) through(Bll) and(B58), the lifting pressure aft of the shock 

? ; I 

'( : in these purely subsonic regions) can be calculated ;-as follows: I 


i/} = ( p u ' p f ° + 4pB j ( 

- + (p 0 B / + g.«] {«,} 


where 



fcol-2 




Expressing Eq.(B76) in terms of downwash via Eq. ( b 5), yields -the desired J 
AIC format given by Eqs.(BlO) and(Bll) : i 


b B } - 5p ba ] {m + b B % v A 


where from Eq.(B5) and(B?6) 


(B-79), 


t“] ■ fjw. “3 • {») I*. “]*F;, J <•-?" 




(B-81) 


Thus, Eqs.(B79‘) and(fs72) display the desired form of the pressure /downwash 1 

. AIC_Ls_w:ith_ prpp.er_c.ons ideration for th e interact ion between th e l ocal i 

subsonic and supersonic regions caused by transmission of disturbances 

through the shock wave and the sonic line. The resultant AIC's are ; 

perturbations of the classical subsonic AIC's (represented by jj^pj - plus! 

coupling terms (represented by P ) ! plus the direct supersonic surface 1 

|7» Ps-A-u — | *- -* | i 

downwash effect ( - k£^ I ) . i 
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APPENDIX C 


PARTICULAR INTEGRAL REQUIRED BY THE LAPLACE TRANSFORM 
SOLUTION FOR THE MULTIlLyERED SUPERSONIC REGION 


I 


I 

In Appendix A, it is shown that in each supersonic layer, the velocity | 
potential obeys the following non-homogeneous ordinary differential equation 
in z after performing a Laplace transformation with respect to X. In layers 


At \ ^ / p f'y* ( v \ I 

0 P zz (s,z) - kI 0 p (k) (s,z) = G k (s,Z) 




where 


4 = {B kS 2 + [2M k Wa k) )2]} 




b -> 


B k s vq - mJ) - , 


if M k < «1 


~ 1 ’ f i<- ) ^ / w ) . 

The fun ct ions <G k (s,Z) arise- from 'initial .'values, of 0^ T ,/ and 0 V ' (on the. 

line ) . 


G„(s,Z) = (S4^> + 


>- SL ' X ’SL 


l 

I 

l J 

I 
I 
I 

(C-2) 

I 

<^3 b 




sonic 


i0. 


00 

[ ’SL 


) + 2i( u ,/a k ) M k 0^ } } "* (C f 5\ 


(k) (k) 

It will now be assumed that the initial values 0 qT and0 x „ are linear in Z 

i on >bL i 

for a narrow layer { ^ Z^\. f | and can be interpolated from numeri r 

cal values at discrete points {& k :k=l,2, --^N }. These discrete numerical J 
values are obtained from subsonic theory as indicated in Appendix B. 
Therefore, f° r ( z k _i /Z Z k ), let 


G k (s,z)' = G^(s) + G k (s) z 


(C-6) 
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where, 

2 = 

Then, Eqs (Cl, C5 through (C7)imply 

/ 


/ 


(G-7) 


«ko(‘) = - t B k <«>sL k) + *i% ) + 


G k (?■) = - {B, ( s 0 J. 


;0 

00 . . (k) 

k V ^9L X 


k - 0 

,(k)i 


(Q-8j 
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) + ■ (< f 9 ] 


and where the linear approximations to' 

I 


fi k L ma 4tsL are 


M (k) (7) - ii\ ^k) , M (k) 

0 SL ^ ” 0 SL- ^SL, 

O 1 


M (k) (7) - _(k) 

^x, SL Z 0 x,SL n 


M (k) z 

^XjSL-l 


= Cu (k) + a (k) 

•' SL 'o 


(G-1Q) 


' (C-lla) 


(e-llb )„ 


(i) • I 

Note that .0l T can be set = 0 on the airfoil surface at the first layer; | 

oi-Jv 1 j 

however, finite values exist for the other layers for transient flow. Next, 
consider solutions to Eq(C-l) ,with(C-6)substituted for the right hand side; '[ 

. i ; 

0p k ^ '(s,z) - K^0 p (;s : ,z) = G k (s) + z \ (s) ' ;(C-12) 

ZZ ' O ■ 1 | 

I 

The general solution of Eq(C-l'2)(or C-6) obtained by the method of variation^ 

j 1 

of parameters, is (ince. Ref 4'3-> p.123) 

! 

0^, k) (s,z) = G k (s,c)dC : 

0 . r~*“\ 

^W( s ih r V^ho G k (=,<)<3c 


^0 


A (0 H j 0 H ) 


(C-13) 


where the homogeneous functions are def ined as solutions of Eq (C-12)' when 
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(s,z) = cosh (K^z) 
(s,z) = sinh (K^z) 


(e-i'4) 

i 

I 

(e-i'5) 


i <b„ ) is the Wronskian determinant: 
H 2 

•~(k) ~(k) 

<t> E 2 

' V ) 

I dz dz 

cosh(K^z) sinh(K k z) 

1 sinh(K^.z) ^cosh^z) 


(C-i6a 


(-Grl'6b 
v . 


, ) = K k (cosh 2 (K k z) - sinh 2 (l^z)) (C-l6c) 


Therefore, for any arbitrary G, (s,z) 


0p( S,z) = — 5 7 - 

P V 


h(Kk z ) 


sinh (Kj^) G R (s,£)d£ 


+sinh(K, z) f ^ „ 

■- I cosh (K^) G k (s, C )d^ (C- 17 ) 

iK u i 


For linear approximation to G k (s,z),j given by Eq(C6), one obtains 

~ k rCosh.(Kz) ~ („\ r -i. 

0 p (s,z) = — p — — G fc £ cosh (I^z) - lj 


+sinh(K k 2 


G k ^ S ^ sinh (Kj^z ) 


iv I 

1 ^ 1 

- |{ cosh (KkZy sinh^z)^ (s)d^} ( 


+ ~Y^ [sinh^ z )J costing )gG (s)dg} (c-jgai) 


Trac'K^wwTr 


FORM 7681-2 




G k (s) < 

Or-.! 


G 4 S > 


2 (- 1 '■ c0S h(K^z)'} - — [:i 1 (s,z) coshCl^z) 

- I 2 (s.,z) sinh(K^z)} (C- 


The integrals, I-^Sjz) and I^^z), are defined as 


V s ’ 2 ) Qsinh(\OK 


I 2 (s,z)l =[ Qcosh (K^cJd^ 
I J 0 


From elementary integral tables 


Ip.Csjz) " ~^2 |^(K k z ) c °sh( K k z)-sinh(K k z)J (C-20a) 


i . I 


Ig.(s , z ). 5.. — _ K^z ) s inh(K k z ) -cos 

\ L 

Substituting (C-20a) and (C-20b) into (C^l8b) yields 

•\ (S) 

0p k \s,z) = — 2 — { -i Hr cosh(K^z)} 

K k 

«(0 ' 


hCK^zl+Uc-J 


K i r 

3— r [cosh(K k z) (i^z) coshCK^z) 

\ ' - L_ ” ( s ^ 

-sinh(K^.z)ji -sinh(K k z) |l-cosh(K k z) 

+(K k z) sinh(K^z)J }' 


(C-21a) 
1 I 


Eq (C21a.) may be further simplified to yield 

~ (s). 


0p k \s,z) = — ^ — [-1 + coshCiyO}^ ij— f(® fc z) 

*k \ ' ' 


21 c 

-sinh(K^,z)j 

k 


(CH’lb) 


& , 

FORM 760 1-2 
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Eq(C21b) represents the function required by the multilayer supersonic 
solution defined in Appendix A. Appendix A also requires and 

evaluated at z = 0 and at z = (Z, - Z, ) = z, . 

x 'k k / k 

first, the derivative of(C21b) is 


~00 

\ (s) \ (s> 

si " h --ir tK k _ 

% 


V)iy 


~ (s) 

°k 

o 


7 (s) 

G k, 


*• 0 p y(s,a) =r#- sinh (Kyz) - — [1 - cosh (K^z)} (C-22b) 


Eq(C22B) represents the interference downwash due to the(particular integral 
in each layer {Z : k = 1, 2:N}. At z = 0 one has from Eq(c21b) 1 


0p k \s,z = .0) s 0 


(C-23) 


At z = 0 from Eq(c22b) implies 


0 p J k) (s,0) = 0 


(C-24) 


At z = z k = (Z fc - Z k _p)Eqs(C21b) and(C22b) yield 

G k^ S ^l 

0p k ^( s 3 z k ; ) = — ~ 2 ~ {- 1 + cosh (K kZk ) } 

!< k 

\ (s) 

- -T-t (\\) - sinh (K^)] 

■Sc I 


(C-25) 


LOCKHtCD 
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(s) 


J k 


I 


(s) 


Z ( k )/, \ "o . , s "1 - . -- 

0 P,z (s ’ Z k) = 


k ; - ~Y — si p h 2“ f 1 ■ cosh (Vk^ (c-26) 

^ I \ , 

I 


The (^quantities defined by Eqs(C23)through (C26)are needed for evaluation of 
layer interface boundary conditions • defined in Appendix A. 


Particular '.Integral for a Subsonic 
External- iFree-'S.tream Disturbance 
I 

.The Laplace transformed disturbance potential is governed by 


0 P (oo) (s,z) - 0^(s,z) = G ( s , z ) 


zz 


where 


£ =_{.B_!sL± 2M «_(-fe)-S— C^--)-}- 


and, in general, 


(C-27) 


-(■C-28) 


= M 1 ; 1 M oo > 1 

CO 00 


= ^ 5 
00 . 


If Moo < 1 it is assumed that initial 


M os < l 1 


(C-29) 


value disturbances are defined on an 


outward extension of the sonic line, and are related to subsonic flow dis- 
turbances in the manner described in Appendix B. If M oo > 1, then no . ! 

initial disturbances need be considered; therefore, the particular integral 
is non-zero only when M » < 1. 

The right hand side of Eq(C27) like Eq(C5) is given by 

i 

l (s.*> ; <bJ(+£> 


(C-30a) 
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- (C-30bi 


It will be assumed that G (s,Z) can be approximated in the form 


G(s;z) =G e _b= ° Z =1G (s,Z )e _b:oZ 

co co . O co IN 


(C-31) 


where for N interior layers 


z = Z - -Z, 


(C-r32) 


Linearly independent homogeneous solutions of Eq(C27) are of the form 


~(«)/ s -Kcoz 

07r '(s,z) = e 

H 1 

-(co), . + K “ Z 

0tj (s,z) = e 


(C-pO 

(C-33b) 


By the method of variation of parameters, one can immediately write the 
particular integral to Eq(C27) in the form prescribed by Eq(C13)* 


0p°°\ s »z) = -e 


+K coZ r I , -K 


(e" KroC ) C?'V) G 6q 


-K 00 z r e (Koo-tto )^' 


G ...dr 


where A is the Wronskian determinant 


I OCKHEED 


FORM 7681-2 
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/ \ r "s K oo Z /% 

••• (t,o )/ 

! Jz 


co - (tfco + K 00 )q 
e d£ 


K co z r ( K co — b co ) r 

' ( G ~ J / e ^ 


/Tv “K co Z 

&“'(G 

<0-K oo | CO jO 


Carrying out the integration yields for (b / K ») 


^(s.c) 




( G ) K co z - (b co + K co )z / \ 

C° ; O r ~s ,e /y-. -K coZ ■ (K oo -b co ) 

*2 K co 1 ,(K co + b co) 'i 6 (e -1 


G - Id co z f ^ —b co z 

°°,o r -e I cte 1 

■2K 00 ( K 00 +b CO ) (K CO -b CO 


G -K oo z .. 

OO ; o e (■ 1 

2K oo l (K »-b co) j 


In the special case b = K co then Eq(C36) yields 


~Go),„ G =°,o r -e" KcoZ ' .ze" KroZ : , 

0 p ( > z ) 2 K oo t 2 K co ! + ^ 


(e 

. K CO -b co . 


/ *G \ -K ooZ 

6 fi (K „z) - i 3 (C- 

2> K„ | 2 ! 

For large K coZ the first term dominates and the disturbance decay rate is ' 

i 

slower than for the case K m ^ b co. Since Eq(C27) also governs the subsonic 
r *** I I 

lifting surface theory (when G ro (s,z) = 0) one could expect the sonic line 

disturbances due to subsonic doublets to behave like e _K Z at large K Jz; I 

therefore, Eq(C39) is probably the relevant solution. ' | i 


ock-hecd 


FORM 7601-2 
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APPENDIX D 


SHOCK WAVE TRANSFER MATRICES 


The following analysis is restricted to normal shock. waves (shock front 
perpendicular to the airfoil surface).. This is. a good approximation to 
many cases of interest in transonic | flow, at least as far as the final 


transition from supersonic to subsonic flow, even when there is a prelim- 
inary supersonic oblique shock compression near the end of the super- ' 
critical zone (above a so-called "supersonic bubble"; see Piercey Reference 
24 Section 11 ..and- Blackwell, - Reference 24 Section 21) f. ‘ , 


The following analysis uses results of Landahl (Reference 5, p. 113) and 
expresses these results in the format required by Appendix B. Landahl I 
writes the relation between velocity potential and streamwise perturbation 
velocity (backwash) across a normal shock as follows: I 


where , 


“ = 2 vlj 


The reduced frequency k^ is defined in terms of a reference length b and 
the upstream flow velocity U. 


k A = f\ 


I " | 

Here we use A^ and to denote conditions upstream and downstream of the 
shock, following Appendix B. The tangential velocity component, W, I 
(upwash) must be continuous across the shock wave; therefore, 1 
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or . 


(— ) 

a*-sH 


./u 

(— ) 

^SH 


■ (D-4b) 

I 


For a single layer, the relations(Dl) to(D4) can be expressed in state I 
vector form (as in Appendix B). Let ^ 


/> 


0 SH 




A\ 


SH 




W, 


u 

SH 


(D-5) r 


where, 


rS 


u„ 


SH 


Then, Eqs(Dl) to(D4) can be written 


a-* 


\ 


X SH } 


(d-6) 


in matrix form 




A 

_ 0 _ 


7 U 

0 


SH 

-■\r 

1 

U SH 

- A u 

w„; T 


1 




SH 


—3a. 


H 


b 


0 


— 1- -0 
-1 0 

0 1 


IL 




/u 

Sw 

- B u 

U SH 

BU 

SH 


> 


(D"7 ) 


w 


The operator g/^x becomes s after Laplace transformation, and"'all harmonic', 
x dependent (barred) quantities become ( ~ ) quantities, dependent on s, 

the Laplace transform complex variable. 


Thus, Eq(D7) transforms to 


- 

- 


f A "] 



- 


— 


r b tt' > 

b 1 

0 


'x U 
0 SH 



-W b 

-1 

0 


'X u 
0 SH 




A tt 







b tt 

s - 1 

0 

< 

~ u 

U SH ' 

► - 


s 

- 1 

0 

1 

— U v 
U ’SH f 




A„ 








0 0 

1 


^ U 
L shJ 

1 

0 

0 

0 


yt u 

SH 


I 


I 


(d-8i) 
i 

i 


I 


Thus, one can write the transfer matrix partitions for change of the state 


:? vbCtor* a'cro'ss' a* normal - sho'c'E a’t the centerline of layer E a s "follows 
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Evaluating the upper left hand partition of (D-ll) yields 


-W b - 1 -W b 1 

s -1 _ s — 1_ 


:-i i 


+ s ) r s -W b 


-ipt/b 1 
s -1 


1 (:W, b+ s) ~- 2 

(_Jai s ) |0 (-s+ia/b) 


Therefore, Eq (D-ll) reduces to 


-• 3 

(k)SH _ 

BA 

L c 


1 -2/(s+ia k /b) o 

o (-s+ia k / b )/( s +ia k / b ) o 


Since the shockwave transfer function is a point function there is no 
coupling .between, layers, in Jihe_larger_matrix_ implied _by_ Eq .(B-26J . 


LOCK HE E D 


FORM 7681-2 
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APPENDIX E 


EXPLICIT SOLUTION FOR TWO INTERIOR SUPERSONIC LAYERS BOUNDED BY AN 
EXTERNAL FREESTREAM AND TECHNIQUE 'FOR GENERALIZATION TO N. INTERIOR IA) 


From Appendix A the Laplace transformed boundary conditions at interface 
and 2 are summarized as follows: At the interface Z = Z^ 


At ;Z = Z 2 : 


(1) 

(2) 

z ( s ’ Z i)/ U i 

= % 

* z 

,(!) 

(2) 

> ( s 5 Z-j^ ) ( 

1= P 

<2) 

(o°) 

z ( s ’ Z 2 } /U 2 

r 

(2) 

(eo) 

' (s,Z 2 ) 

/N>. 

r P 


J 2" » 


(s,Z 2 ) 


Using Eqs (15) through ( 17 ) of appendlix A, Eq.(’E-l) becomesfat Z=Z, ;z=z =Z n )1 

I -Lilli 


[c, sinh(K- L 


1 1 K-D- 

Z l) + D 1 COSh (^ Z l)j + u/pz < B * Z 1> = -U^ 


~P ! (U x s + iw^C^ cosh( KjZj) + D x sinh( JJ +^ 1 ^(s,z 1 )|= 

- -p 2 (U 2 s + i u) C 2 


At the external freestream interface (Z = Z 2 , z = = (Z^Z^)): 


i 

C 2 sinh(K 2 : 


' 2 ' ^ 2 


osh (K g z 2 

p (2) 

} J + U 2 ^Pz ( s ’ z 2 ) = 

n 

8«|- 



So C co + ^P z (s,Z = 0)| 

2 z 2 ) + d 2 

sinh(K 2 z 2 )l +#£ 2 )(s,z g ) 


7 (-) r 

-=--p=°(-uj 

s -+- ±Cd)--|-C oo -+-^- -(s,0)| 


(E'-2c 


LOCKHEED 


FORM 7681-2 
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Next, we write the Laplace transformed airfoil surface normal velocity 
(upwash) boundary condition I 


^(s) 


k a + 


(E-3) 

I 


Eqs.(E2a-E2d)and(E3) provide five equations for the five unknowns, , D ,| 

C 2’ D 2’ C oo* 1 

I 

I 

Each of the interface pressure and streamline slope continuity boundary I 
conditions can be put in the form of a transfer matrix. 


7 «Y + P ^TR). R? #| m 


where for Z = Z., k = j or j + 1 we define vectors 

I 


i F m 

\ p i l 


such that 


MT= 

R i) 


i / u k ^ k) (s,z. ) / a z^] 


Pp k \s,Z.) 

j 


(e-4) 


(B-5) 


P p . = '" p k^ U k S + \ ^ 


The vectors 


m 


(e-6) 


define forcing functions arising from the particular' 


integrals (defined in Appendix C) which, in turn, are caused by the sonic. 


line initial value disturbances. 


The external freestream interface conditions (E-2c) and E-2d) can be written 


in similar format 


Hff + P’Jftj * I 9 '"! c » ft? <E f 7) 

i' ' 1 

where, (o=) implies the use of external freestream steady flow field conditions 
for fluid property (l^, evaluations . The two by two transfer matrices 


Loc tfffrso- 


FOfcM 7681-2 




^A (k) ] and [b«J are defined as follows from Eqs (E-2) 

A U = r sinil (K k\ J 

k 

A i2 J = r oosh (K A> 

k I 

A 21^ = " p k^ U k S + ^ c oshi( K^z^.) 
k 22 = “ P k^ U k + ^ sinh ( K k z k ) 


^ . 0 


B ^ +1) * A 12 +1 | ^ S,Z k+l = 


,(k+l) _ ,(k+l) 


l 21 l< s ’ z k + l “°)=-P k+ l(\ + i s+ i“) 


B<| +1) =Ai k+1 !( S> k wl = 0) = 0 


Thus, the elements of and Fb ^ k fj satisfy 

which is helpful for computation. 1 


a recursion relationship 


The column vector {e^} defined in Et * ( E ~^ is associated with the inter- 
face between the outermost interior layer and the external freestream( valid 
whether is either supersonic, or subsonic). ^(3^*^Jis thus defined Iby * 

Eqs (E-7) and(E-2d) by an obvious generalization from 2 to N supersonic j 
layers. I I 


I 


l 

l 

i 




^li 1 ) 

i bH+bH'’ 

\ 21 22 i 


(E-10) 


The superscript (oo) implies substitution of freestream fluid properties for 
the (N + l) layer in Eq(E9) • Thus, I 


B n ) " 0 > $2 = - k =A 


B 21^ = " P“( U I S + ^ 5 B p2^ = 0 


(E-ll) 


Therefore, (E-lOijand(E-ll) yield 


-K /U 
00 ' 00 


— cl 1 - t" p “fe s - + ^f- 

The two v by xfwo transfer matrices and P +1) ] 


(E-12) 


The two^by /two transfer matrices LT'J and 

jlJ J allow one to easily 

solve for the coefficients -j^j- including especially and D-j^ which define 
the surface pressure. First,. it is necessary to eliminate C by use of Eg 

| 00 I ^ 

(E-7) which can be re-expressed as a constraint between C -and D (or C I 

N j 

and D in the general case of N interior layers). . 


C N . a M D N + C NP 


where, from Eqs (E-2c) and E-2d) 


(E-!3) 


A 2 } C N + A 12 ) V (1/U N ) '^ S ’ Z N ) ^ + 

00 I 

A 21 C N + A 22 " p N^ U N S + ^ ®P^ S ? Z N^ = "P 00 ^® 3 + \ s >0)) ! (E-l4b) 


! . O C KHifCtS 


FORM 7681-2 





Solve Eqs (E-l4) for 


a .(N) + A (N ]d + p (N) (s z ) 

A 21 C N A 22 D N P P,N ’ V 


= Pp.N ( s >°) - Pco'(U ro s + ito) 

-A (N) C - A ( ^D ~ 
A 11 L N 12 Tff 


^ r 

k . L 


fA (») ,(nK 

' P,N ' P,N 


(E i 15) , 


Therefore, (E-15) can be written 


(N)„ (N) ^ (») 

Yu C N +Y 12 , d n •'% 


(E-'l6) 


where, 


& - 4? - *<l- ♦ ^ .‘if 


(E 1 17a) 


(N) _ a (n) ,1 „ . v j» M 

Yl2 A 22 " P ^Y» S K A 12 


,(“) - ^ (“) 


(E-S17D) 


-Vi(-) -*(N)w4 


«r - w>°> - + r<*& f) 


and the particular integral contributions to interface streamline slopes are 
defined as I ' 


*p|k- = uj^ ^ S >V 

. (») I 1 — (<=) 

^ TD TVT “• TT ^ T>r7 ' S 5^ ) 


P,N U ^Pz' 
• 00 


(E-19) 


From Eqs (E-13) and(E-l6) through(E-l^) 


(IN)/ (N) 

“nn ■- Y i2 ' Y 11 


( E— 20) 


’ N,’P % “ Y-Yn- 


”(-E— 2l) 

; I 
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Then, 


3 2 = Tl,2)! < D 1 + AD ljP ) 


(e-28) 


where &B is defined by 

if 


-a (1) 


Agg P>1 • . 


2,P .' 


( E- 29) 


From Eq ($-26) 


(1,2) _ (1) (2) J (1) (2) 

E 1 - I A 11 B 12 + \ *22\l 

(1,2) _ (1) (2) (IV (£ 

E 2 I Ag 1 B 12 J A 22 °^22 B 21 


(E-30a) 


(E-30b,) 


‘ E i 1,2) - ^21° 


(E-30c)^ 

I 

(E-r30i) 


Next, solve the first' row of Eq (E-2U) for C , since is already known I 
from Eq (E-3) 1 


f — -i- i + t 

C ! " E l °2 + AE 1 C 2,F + X A 




,(1,2) 


|pet (d i + 4c i,p ) ♦^’S.p ♦ f 3£) 


Some general simplifications occur as 'a result of results in Appendix C | 
concerning the particular integrals ^^^(s,z). These are expressed as I 

l-L . 


T pfi <s ’V = 8 p>’°> =X ? p' ,( s ’ Q) =0 


^ 2) (s,0;> =^ 2) (s,Z 1 ) -= 0 


( E-3 3 ) 

i 

( E-B^ ) 


Pp 2 ^( s > z ^= Pp^( s >0) = -P 2 ( U 2 S + i^) ^p 2 {(s,0) = 0 ( E-35 ) 



FORM 7681-2 




Likewise, Eq(E-3) for simplifies ( 


e^Cs.z - z = o = ^-?^[s,o) = o 


(e-'b6) 


Pp 1 ^(s,0) = Pp|i^(- s »0) = -P 1 ( U 1 S + i w)0p 1 \s,O)=o;(E-37) 


Then from Eqs (E-3) and (E-36) 


D 1 = \( s V K i = \ 


(e-38) 


which agrees with Reference 2, p. 366, for the Laplace transform solution 
to the Garrick and Rubinow case of an unbounded uniform supersonic stream. 


Then using Eqs (E-34‘) to (E-38) one obtains from (E-29) that aD-j 
simplifies to 

-ad— - I (1) 6- (1 > +1 : (1) D (1) + 

“i,p - \ r 6 P,i + 1 a z2 P P,1 + “2 C 2,P 

p (l,2) 

C 1 - JIT2T < D 1 + * D 1,P> + '-4 1,2)c 2,p 

t- 2 I 

_ (T (1) (1) + (1)5 (1 \ 

U A 11 6p ,l J A lT V P,l } 


■(E-39)”f 


(e-4o) 


In the 2 interior layer case, N = 2 and Eq (E-21) yields 


C = 1 

C 2R . S 2 /Yu 


where from (E-l8)' 


(e-4i) 




«p?2 = ?k )(s>S 2 5 


(e-43) 




(e-44) 





(E-j45) 


(e-46; 

I 

I 

(e 7 47) 

i 


i 

i 

«» 


(e-48) 



:(Er50) 

I 


'(E-51a) 

l 

(E-518 ) 

I 

I 

(E-51c) 

(Et51<0 



1 
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~\ Forming the matrix one finds ( 

-U 1 sinh(jK 1 z 1 ) -coshCK^) 


( 1 ) 




K, 


U, 


P l( U i s + ^ 

sinl^K^z ) 


(E-52aiO. 


<j-) oosh(K 1 z i ) 

'2 1 (s)sinh(K_ L z 1 ) -coshtl^z^) 
-*2 1 (s)cosh(K 1 z ) sinh(K z ) 


p,(U s+ioj) 


where we define the Laplace transformed interface impedance in layer 1 asj 


(E-52b ) 


I 


'Ms) =; -P 1 (U 1 S + ito)U n /K 


'1' 1 


(e- 53 )l 


Airfoil Surface (Pressure in the 
Supercritical Region 


As defined in Appendix B the local surface pressure in Ay(where > l) 


A. 


~ M ~ A Uv n s 
P a ° = , P -(Z = 0) 


u- 


is 


-P].^! 8 + ic ° ^ (C ]L coJ sh (K 1 z) + DjSinhCl^z) +^ 1 ^(s,z).^ _ 0 


- + i|o) (C x +?p 1 \s,0)) 


(E-54'): 


But from Eq (E-ST) 0p^(s,O) = 0; therefore. 


A A i 

Pg U = P U (z = 0) = -p 1 ((U 1 s + ico) C-j 


(e-55 ) 


i 

Using Eqs (e40 -E53) for C, yields 

J. i 
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E. 


-(- I72 -)- 


V(s) 


( 1 , 2 ) 


C = (12) ( ~K ) + - 1 ,' rtD . + AE '(l,2) r 

C 1 E^ ,2 > K 1 E | (1,2) AD 1,P AE 1 C 2,P 


(u, s ' V i a) <-' 2 1 ( s )f i,,h (K 1 z 1 ) 5 PA + «E56) 


( 1 ) 


p l'~l 

Now the ratio /]£ is the 

From Eqs (E-30a,Er30b,E-52a. & E-52b 

(1,2) JL 


I 


E_, 


quantity called E (s) in Appendix A. 

^ X (l 2 ) ( 12 ) 

,we first define E^ ' 5 ' and Eg' ’ / . 


p 1 (U 1 s X T ico) ^(Zi(s)Bj 2) sinh( K iz) - O^r^'oshCVi^fE-ST) 

1 1 I 

pUuA to) ^ 1 ( s Ph2 2)cosh (h z 1 >“22 B 2 [^ i “ h (Vl»M E -58) 


( 2 ). 


( 1 , 2 ) 


’ l'~l 


/ 

From Eqs (E-8, E-9,' E-45-E-51') 


vif = ~P p(Up s ” + -ico)cosh(K p z ? ) - p a> (U B S + ini) sinh^Zg) {E- 59 ) 


V K 2 


,(■ 2 ) 


t 


K U_ 
00 2 


^ K 2 


Y nP = -pp(U ? s- + ito)sinh(K p z p ) -^(U^s + ico) p-p- coshfKgZg) ' (E-60) 


’12 K 2 2 


K U 0 

03 2 


4- » 




‘ ! 


a 


22 


-v < 2 > y v ( 2 ) 
^12 ( Y n 


'(sinh(K 2 z 2 ) + p 


p /U s + ico) U K, 

’ nn I rr\ CO , 


'03 “2 


2 V + fcyyij 


-cos 


h ( K 2 z 2 V ) 


• 0 ; (U.'s +‘ ico)U K 0 7 ■ " 

(cosh(KgZ 2 ) + p 2 (Ug S + ico)U p ir Sirlh ^ K 2 Z 2^ 


J 2 » 


- <sinh(K 2 z 2 ) + ^(s^coshCKgZg) j 



ci 


%pj 


22 


(cosh(KgZ 2 )+ ^co^j sinhCKgZg) j 


(E- 62) 
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(sinh(K 2 z 2 ) + 2 (s) cosh(K 2 z 2 ) 

22 _ (cosh(K 2 z 2 ) +'Z co 2 (s) sinh(K 0 z 0 ) 


(1E-63) 


2 2 ‘ 


where the impedances, ^(s) and *2 (is) are defined by analogy to Eq (E-57) 

d 00 1 ^ 


(s) 


= -P 0 ('U 0 s + i^)U_/K. 


and the ratios 


= -p (IU s + iaj)U /K 


■2 2 j 1 (8) =^,(31/2, (s) 


From Eqs (E-57 - E-66) and Eq (E-9) 


( 1 , 2 ) 


( 1 , 2 ) 


P i (U 1 s + 1 m) Ci l ( ' S) uf Slnh(K X Z l ) * «22 P 2 (U 2 S + to, )c° sl '(K 1 q)'!l 


(E-64) 

I . 

Ce- 65) 


(e, 66) 


(e-67) 


K_ 


(1E-68) 


V"*' * p j tqiy M COSh(K 1 z 1 )' a 2f2( U 2 s + HoJsinhfK^zp ) 


We now form the ratio 


( 1 , 2 ) 


e 2 i ( 3) \( s) 


, (1,2) 

7e 2 (s) 


^ 2i ( s ) = sinh(K 1 z 1 ) . o; 22 cosh( )) 

("\/ 2 2^ COSh ( K l Z i) +a 22 Sinh (^L Z l)^ 


(IE- 69) 
I 


= _-(tanh(K 1 z 1 ) “ \ f x ( s )» 22 ) 

(i - o 22 'Z 2 ^ 1 (s)tanh(K 1 z i )) ; 


(E-70) 


t I 
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Now, from Eq (E-63 ) 5 0! 22 can wr iti en in the form 

-( Z a>,2 (s ) + | taah (K 2 z 2 )) ^ / N 

a 22 (1 +'2 oo ^tdnh(K 2 z 2 )).' V Sj 

where H (s) is the quantity defined in Eq(A-32)Appendix A. 

If one defines as in Eq (A-29) of Appendix A 


H 21 (s) 


Then Eqs (E-70 - E-72) yield 




\l(s) - . 


(tanh(K 1 z 1 ) - < B gl (s)) / 


(E-j7l) 


(l - T^ 1 (s)tanh(K 1 z 1 jj' 

This is precisely the same as Eqs (^27 ) through(A~33) of Appendix A. 
To further establish the correspondence with Appendix A we note that 
(A-33) of Appendix A equals Eq (A-71) and from Eqs (A-27) and (A-30) 
Appendix A and Eq (A-78) it is seen that 


( E 7 2) 


(E-j73 ). 


Eq 

of 


K 2 i<s) - 


(fS 


21 (s) - tanh(K 1 z 1 ) / ) 


D 2 i( s ) = ,( 1 - 'H 21 (s)tanh(K 1 z 1 ) ) 


* 


(e-'7M 


(E-,75) 


where from Eqs (E-66 - E-71? E-72 and E-73) 
H 2 l(s) - 


P 2 (U 2 s + iqj j s ) + tanh(K 2 z 2 ) ) 


P 1 (U 1 s + iio) U-jK^l + 'Z c05 2^ S ^ tanh ^ K 2 Z 2^ 


(e-76) 


Particular Integral Contribution to 
Local Supersonic Surface Pressure 


It can be seen from Eqs (E-54 - E-56 E -69 - E- 76 ) that the term 

T5 21 (s) W a (s)/K^ represents the direct contribution of the local surface! 
upwash, y w a ( s ) in the local supersonic region Ay The remaining terms in 


, i | 

Eq (E-56) represent surface pressure effects due to sonic line disturbances. 

Following Eq (B-42) of Appendix B. j j 

93 -> 


'CKHEED 






V 

P (z = 


' I 

0) } ' { P ^ U } = twA AU ] : { M A + iW A P ) + f P,S U ) < E "77> 


Eq (E-77) represents, in matrix form the numerical evaluation of the conyolj- 

ution integral, which is the required inverse Laplace transform of Eq (E-55 ; )$ 

| 

with Eq (E-56) substituted for C p . I 1 


Eq (E- 56 ) can be split into two terms 


c i = c i,h : + c i,p 


(e-78) 


where the first term, is the homogeneous solution for the multilayer super- 


sonic region. 


C 1,H E 21p^ W k( s )/ K l ( 

and the second term of (E-78) is the^ particular integral contribution. 

c i )P .= e 21 (s) 4 d 1;P+4 Ei (1 > 2) c 2! p 


(E-79) 


P 1 (U ] _s + ito) 


(-Z^s) sinh(K 1 z 1 ) 6 p ^ + p p ^ C osh(K 1 z 1 )) (E- 80 ) 


1 1 

By comparing Eq (A-54) with Eq (A-22)pf Appendix A and Eq (E-77) above, it 

can be seen that the Laplace transform of the particular integral contribu- 
tion to the surface pressure in the Ideal supersonic region is I 


APp,i( s j 0 ) - 1 L (p p ) 


(E- 8 l)j 


=! -p^UpS + i^)C ljP 


Since, by Appendix C, there is no direct contribution to the airfoil 1 

1 | 

surface downwash from the particular integral in the first layer, we have; 

! 

J&pi 1 1 ( s .> 0)_ ^AW Ap (s) _=_ 0 _(E r 82) H 
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-Then- from -Eqs ( E-55 , • E- 78 " and E-79) \ 

A.. I V (s) 

P U (s,z = 0 ) = -^(UjS + iai)(E 21 (s) -|^4 + C 1 ^ p ) 


(E-83) 


Therefore, by the convolution principle 

• x 


I 


P U (x;,z = 0) = / E« (x - |);Ap^(?)d? + L" x (-p (U s + ito)C p ) (E- 84 ) 


GR 


1 V 


21 


•v 


l'~l 


1 ,P' 


A p GR (f) = L _1 (-pl(U lS + itojw^s)/^) 


(E- 85 ) 


An alternative form of " Eq (E-84) shows the direct comparison with Eq (E- 83 ) 


x 


U 


Vu 


p u (x,z =0) =/ G m u (x- ? ) W a ( 5 )d? + L" ( -p -j ( U-, s + i 40 )C n p ) (E- 87 ) 


'0 


l'~l 


1 ,P 


where by the Laplace transform inversion theorem 

h A „ 1 r c+1 " 1 .fv.i 


2 rFi 


.( -p r l 


c - 1 °° 


Thus, it is seen that Eq (E-77) is a 


• U I s_+_ixo). ~^= )-e_ v “*'ds (,E= 88 ). 


matrix statement of a numerical inte- 


gration scheme which is equivalent to the analytical solution expressed liy 

I 1 

Eq (E- 87 ).. Therefore, Eq (E- 87 ) ' represents the basis for calculating the J 

» 1 I 

matrix elements required by Eq (B-4l) of appendix B. I j 

! 1 

! 1 I 

We next consider the inversion of C j p . First, Eq (E- 80 ) must be simplified. 

* -L 5 ir • 

Using Eqs (E--52a,E-52b,E-51 and E-29) one finds 

(1) (1) (1)|~ (1) (1,2) - -i 

AD 1,P X A21 6 P,1 J A22 j P P,l AE 2 C 2,P 


Pl (P lS 4 to j (-^(^cosM^p »p^+ : | 

: ' 1 1 

: (E- 89 )! 


(1) 


+ AF 

+ ae 2 c 
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From Eqs (E-4l - E-59) 

1 

1 

1 

C 2,P ■ g 2 < “ )/ Vli 2> ■ 

(E-90) 

(p p j o) ( S ,o) -~P p , 

( 2) (s,z 2 )) ! 


-p 2 (U 2 s + ito)(cosh(K 2 z 2 ) + Z^^sinh (K^) ) 
From Eqs (E^^b, E-56b and E-9) 


(1,2) _ (1) (21) 

4E 1 " Z A12 B 21 


P 2 ^2 


(U 0 s + ioo) 


— 1 — — ... \ (-cosh(K n z )) 

P l( U l S + lt0 ) 1 1 


(E-91) 


i Substituting Eqs (E-89 j E-90, and E-9I into E-80) provides the particular 
| integral contribution to the local surface pressure in the supersonic regibh 
. which is represented by Eq (A-20) ofl Appendix A 


£P Pjl ( s > 0 )'' = "P 1 ( u 1 s + ic °) C i 5 p 

= -E 2 l^ s ) cosh(K 1 z 1 )e p ^+ sinh(K 1 z 1 )p^(s,z 1 ) ) 


(E-92) 


_ (_^ 1 posh( Ki z 1 )ft p ^)+^ ? ^)(s,z i )Binh(K 1 z 1 ))+(-E 21 (s) S inh(k 1 zj L ) + 

((Pp^ -k^) -'2 s («p < 2 ) - e < 2) )) 

+ cosh(K 1 z 1 ))' ^ co ’ h( K z i )’+'2 0 sinh(K 0 z 0 )j ’ 


i I 

(e-93)! 


2 2 “<*,2 - ^2 2 ' 

We next simplify Eq (E-93)- It is convenient to represent Eq (E-93) as 
follows .* ! 

(1) . „ (1) 


4p Pj i( s 5 °) = A 


rp,i + b i p p,i * 




(e-94 )| 

i i 

■ i ! 

.. j 
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where 'by comparison of Eqs (E-93 and E-9*0 

\ = Z^s) (1 + E 21 (s)) cosh (K 1 z 1 ) (E-95&) 

\ = - sinh(K 1 z 1 ) (1 +E 21 (s)) (E-95b) 


_ (cosh(K 1 z. 1 ) - E 21 (s)sinh(K 1 z 1 )) 
A 2 - (cosh(K 2 z 2 ) + Z^ 2 sinh(K 2 z 2 )) 


(E.-95c) 


Eqs (E-92) to (E-95) display the effects of relfections , back to the airfoil 
surface pressure, of the contributions of the particular integrals in 
regions (l, 2 and «) to the streamline slopes and pressures at the layer 
interfaces at z = z^ = Z^ and at Z = Z (where z 2 = Z 2 - Z^). The last 
term in Eq (E-9*0 displays the impedance mismatch feature encountered in 
crossing the outer layer adjacent to the freestream flow. The inversion of 
the Laplace transforms implied by Eqs (E-93) to (E-95) will be discussed in 
a later study.,.'. 
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